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ABSTRACT 
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We report measurements of the mass density, 51m, and cosmological-constant energy density, of the uni- 
verse based on the analysis of 42 Type la supernovae discovered by the Supernova Cosmology Project. The 
magnitude-redshift data for these supernovae, at redshifts between 0.18 and 0.83, are fit jointly with a set of su- 
pernovae from the Calan/Tololo Supernova Survey, at redshifts below 0.1, to yield values for the cosmological 
parameters. All supernova peak magnitudes are standardized using a SN la lightcurve width-luminosity relation. 
The measurement yields a joint probability distribution of the cosmological parameters that is approximated by 
the relation 0.8 17m -0.6 fl\ ~ -0.2 ±0.1 in the region of interest (f^M <, 1.5). For a fiat (J7m + ^a = 1) cos- 
mology we find Sljjj" = 0.28^Qg (ler statistical) ^{jl (identified systematics). The data are strongly inconsistent 
with a A = fiat cosmology, the simplest inflationary universe model. An open, A = cosmology also does not 
fit the data well: the data indicate that the cosmological constant is non-zero and positive, with a confidence of 
P(A > 0) = 99%, including the identified systematic uncertainties. The best-fit age of the universe relative to the 
Hubble time is $ at = 14.9±J; 4 (0.63/7i) Gyr for a flat cosmology. The size of our sample allows us to perform a 
variety of statistical tests to check for possible systematic errors and biases. We find no significant differences in 
either the host reddening distribution or Malmquist bias between the low-redshift Calan/Tololo sample and our 
high-redshift sample. Excluding those few supernovae which are outliers in color excess or fit residual does not 
significantly change the results. The conclusions are also robust whether or not a width-luminosity relation is used 
to standardize the supernova peak magnitudes. We discuss, and constrain where possible, hypothetical alternatives 
to a cosmological constant. 
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1. INTRODUCTION 

Since the earliest studies of supernovae, it has been suggested 
that these luminous events might be used as standard candles 
for cosmological measurements (Baade 1938). At closer dis- 
tances they could be used to measure the Hubble constant, if an 
absolute distance scale or magnitude scale could be established, 
while at hi gher redshifts they could determ ine the deceleration 
parameter (Tammann 1979; Colgate 1979). The Hubble con- 
stant measurement became a realistic possibility in the 1980's, 
when the more homogene ous subclass o f Type la supernovae 
(SNe la) was identified (see Branch 1998 ). Attempts to measure 
the deceleration parameter, however, were stymied for lack of 
high-red shift supernovae. Even after a n impressive multi-year 
effort by N0rgaard-Nielsen et al. (198"^ ), it was only possible to 
follow one SN la, at z = 0.31, discovered 18 days past its peak 
brightness. 

The Supernova Cosmology Project was started in 1988 to 
address this problem. The primary goal of the project is the 
determination of the cosmological parameters of the universe 
using the magnitude-redshift relation of Type la supernovae. In 
particular, Goobar & Perlmutter (1995) showed the possibil- 
ity of separating the relative contributions of the mass density, 
£1m, and the cosmological constant, A, to changes in the expan- 
sion rate by studying supernovae at a range of redshifts. The 
Project developed techniques, including instrumentation, anal- 
ysis, and observing strategies, that ma ke it possible to systemat - 
ically study high-redshift supernovae ( perlmutter et al. 1995a ). 
As of March 1998, more than 75 Type la supernovae at red- 
shifts z = 0.18-0.86 have been discovered and studied by the 



1996 



Supern ova Cosmol ogy Project (Perlmutter et al. 1995b 
1997g@gg[l9983) 

A first presentation of analysis techniques, identification of 
possible sources of statistical and systematic errors, and first 
results based on seven of these superno vae at redshifts z ~ 0.4 
were given in Perlmutter et al. (1997e; hereafter referred to as 
"P97"). These first results yielded a confidence region that was 
suggestiv e of a flat, A = universe , but with a large range of un- 
certainty. |Perlmutter et al. (1998b| ) added a z = 0.83 SN la to this 
sample, with observations from the Hubble Space Telescope 
and Keck 10-meter telescope, providing the first demonstration 
of the method of separating £1m and A contributions. This anal- 
ysis offered preliminary evidence for a low-mass-density uni- 
verse with a best-fit v alue of = 0-2 ± 0.4, assuming A = 0. 
Independent work by |Garnavich et al. (1998a ), based on three 
supernovae at z ~ 0.5 and one at z = 0.97, also suggested a low 
mass density, with best-fit i l^ = -0. 1 ± 0.5 for A = 0. 

Ferlmutter et al. (Lyytsc) presented a preliminary analysis 
of 33 additional high-redshift supernovae, which gave a confi- 
dence region indicating an accelerating universe, and barely in- 
clu ding a low-mass A = cosmology. Recent independent work 
of Riess et al. (1998), based on 10 high-redshift supernovae 
added to the Garnavich et al. set, reached the same conclusion. 



Here we report on the complete analysis of 42 supernovae from 
the Supernova Cosmology Project, including the reanalysis of 
our previously reported supernovae with improved calibration 
data and improved photometric and spectroscopic SN la tem- 
plates. 

2. BASIC DATA AND PROCEDURES 

The new supernovae in this sample of 42 were all discov- 
ered while still brightening, using the Cerro Tololo 4-meter 
telescope with the 2048 2 -pixel prime-focus CCD camera or the 



4 x 2048 2 -pixel Big Throughput Camera (B ernstein & Tyson 
1998). The supernovae were followed with photometry over 
the peak of their lightcurves, and approximately two-to-three 
months further (^40-60 days restframe) using the CTIO 4-m, 
WIYN 3.6-m, ESO 3.6-m, INT 2.5-m, and WHT 4.2-m tele- 
scopes. (SN 1997ap and other 1998 supernovae have also been 
followed with HST photometry.) The supernova redshifts and 
spectral identifications were obtained using the Keck I and II 
10-m telescopes with LRIS ( |Oke et al. 199% and the ESO 3.6- 
m telescope. The photometry coverage was most complete in 
Kron-Cousins /?-band, with Kron-Cousins /-band photometry 
coverage ranging from two or three points near peak to rela- 
tively complete coverage paralleling the /?-band observations. 

Almost all of the new supernovae were observed spectro- 
scopically. The confidence of the Type la classifications based 
on these spectra taken together with the observed lightcurves, 
ranged from "definite" (when Si II features were visible) to 
"likely" (when the features were consistent with Type la, and 
inconsistent with most other types). The lower confidence iden- 
tifications were primarily due to host-galaxy contamination of 
the spectra. Fewer than 10% of the original sample of super- 
nova candidates from which these SNe la were selected were 
confirmed to be non-Type la, i.e., being active galactic nu- 
clei or belonging to another SN subclass; almost all of these 
non-SNe la could also have been identified by their lightcurves 
and/or position far from the SN la Hubble line. Whenever pos- 
sible, the redshifts were measured from the narrow host-galaxy 
lines, rather than the broader super nova lines. The li gh tcurves 
and several spectra are shown in Perlmutter et al. (1997e, 



1998c, 1998_bJ); complete catalogs and detailed discussions of 



the photometry and spectroscopy for these supernovae will be 
presented in forthcoming papers. 

The photometric reduction and the analyses of the 
lightcurves followed the procedures described in P97. The su- 
pernovae were observed with the Kron-Cousins filter that best 
matched the restframe B and V filters at the supernova's red- 
shift, and any remaining mismatch of wavelength coverage was 
corrected by calculating the expected photometric difference — 
the "cross-filter .^-correction" — using tem plate SN la spectra, 



as in Kim, Goobar, & Perlmutter (1996|). We have n ow re 
calculated these K corrections (see [Nugent et al. 1998 ), using 
improved template spectra, based on an extensive database of 
low-redshift SN la s pectra recently made available from the 
Calan/Tololo survey ( Phillips et al. 1998). Where available , 



IUE and HST spectr a (Cappellaro, Turatto, & Fernley 1995 



Kirshner et al. 1993) were also added to the SN la spectra, 



1990N, 1991T, 1992A, and 1994D in: Kirshner & Oke 1975 


Branch et al. 1983 


Phillips et al. 1987; Jeffery et al. 1992 


Meikle et al. 1996; 


Patat et al. 1996). In 


Nugent et al. (1998) 



we show that the /^-corrections can be calculated accurately for 
a given day on the supernova lightcurve, and for a given super- 
nova lightcurve width, from the color of the supernova on that 
day. (Such a calculation of K correction based on supernova 
color will also automatically account for any modification of 
the K correction due to reddening of the supernova; see N ugent 
et al. 1998| . In the case of insignificant reddening the SN la 
template color curves can be used.) We find that these calcu- 
lations are robust to mis-estimations of the lightcurve width 
or day on the lightcurve, giving results correct to within 0.01 
mag for lightcurve width errors of ±25% or lightcurve phase 
errors of ±5 days even at redshifts where filter matching is 
the worst. Given small additional uncertainties in the colors of 
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supernovae, we take an overall systematic uncertainty of 0.02 
magnitudes for the K correction. 

The improved Zf-corrections have been recalculated for all 
the supernovae used in this paper, including those previously 
analyzed and published. Several of the low-redshift supernovae 
from the Calan/Tololo survey have relatively large changes 
(as much as 0.1 magnitudes) at times in their ^-corrected 
lightcurves. (These and other low-redshift supernovae with 
new /f-corrections are used by several independent groups in 
constructing SN la lightcurve templates, so the templates must 
be updated accordingly.) The -corrections for several of the 
high-redshift supernovae analyzed in P97 have also changed by 
small amounts at the lightcurve peak {AK(t = 0) <; 0.02 mag) 
and somewhat larger amounts by 20 days past peak (AK(t = 
20) <; 0.1 mag); this primarily affects the measurement of the 
restframe lightcurve width. These /^-correction changes bal- 
ance out among the P97 supernovae, so the final results for 
these supernovae do not change significantly. (As we discuss 
below, however, the much larger current dataset does affect the 
interpretation of these results.). 

As in P97, the peak magnitudes have been corrected for the 
lightcurve width-luminosity relation of SNe la: 



m B on = m B + A con (s), 



(1) 



where the correction term A corr is a simple monotonic function 
of the "stretch factor," s, that stretches or contracts the time axis 
of a template SN la lightcurve to best fit the observed lightcurve 
for e ach supernova (see Perlmutter et al. 1995a, 1997e; K im 
et al. 119981 



upernova (see i-'erlmutter et al. iwaa, iw/e; K im 
t poldhaberet alj|1998fc and cf. |Phillips||1993|; R iess, 



Press, & Kirshner 1995, |1996| T A similar relation corrects the 
V band lightcurve, with the same stretch factor in both bands. 
For the supernovae discussed in this paper, the template must 
be time-dilated by a factor 1+z before fitting to the observed 
lightcurves to accoun t for the cosmological lengthening of the 
super nova timescale dGoldhaber et al. 1995|; Leibundgut et al. 
1996a; |Riess et al. 1997a[ ). P97 calculated A con (s) by translat- 
ing from s to Amis (both describing the timescale of the super- 
nova event) and then usi ng the relation betwee n Am 15 and lumi- 
nosity as determined by Hamuy et al. (1995). The lightcurves 
of the Calan/Tololo supernovae have since been published, and 
we have directly fit each lightcurve with the stretched template 
method to determine its stretch factor s. In this paper, for the 
light-curve width-luminosity relation, we therefore directly use 
the functional form 



A COIT (s) = a(s-l) 



(2) 



and determine a simultaneously with our determination of the 
cosmological parameters. With this functional form, the super- 
nova peak apparent magnitudes are thus all "corrected" as they 
would appear if the supernovae had the lightcurve width of the 
template, 5=1. 

We use analysis procedures that are designed to be as sim- 
ilar as possible for the low- and high-redshift datasets. Occa- 
sionally, this requires not using all of the data available at low 
redshift, when the corresponding data are not accessible at high 
redshift. For example, the low-redshift supernova lightcurves 
can often be followed with photometry for many months with 
high signal-to-noise ratios, whereas the high-redshift supernova 
observations are generally only practical for approximately 60 
restframe days past maximum light. This period is also the 
phase of the low-redshift SN la lightcurves that is fit best by 
the stretched-template method, and best predicts the luminosity 



of the supernova at maximum. We therefore fit only this period 
for the lightcurves of the low-redshift supernovae. Similarly, at 
high redshift the restframe B-band photometry is usually much 
more densely sampled in time than the restframe V-band data, 
so we use the stretch factor that best fits the restframe B band 
data for both low- and high-redshift supernovae, even though at 
low-redshift the V-band photometry is equally well sampled. 

Each supernova peak magnitude was also corrected for 
Galactic extinction, Ar, usi ng the extinction law of C ardelli, 
Clayton, & Mathis (1989), first using the color excess, 
EXjj-VOs F&p, at the supernova's Galactic c oordinates pro- 
vided by |Schlegel, Finkbeiner, & Davis (1998 ) and then — for 
comparison — u sing the E(B-V)b&h value provided by B urstein 
& Heiles (1982, |1998| ). A R was calculated from E(B-V) us- 
ing a value of the total-to-selective extinction ratio, TZr = 
Ar/E(B-V), specific to each supernova. These were calculated 
using the appropriate redshifted supernova spectrum as it would 
appear through an 7?-band filter. These values of TZr range from 
2.56 at z = to 4.88 at z = 0.83. The observed supernova colors 
were similarly corrected for Galactic extinction. Any extinc- 
tion in the supernova's host galaxy, or between galaxies, was 
not corrected for at this stage, but will be analyzed separately 
in Section 4. 

All the same corrections for width-luminosity relation, K 
corrections, and extinction (but using TZb = 4.14) were applied 
to the photometry of 18 low-red shift SNe la (z < 0.1 ) from the 
Calan/Tololo supernova survey ( |Hamuy et al. 1996 ) that were 
discovered earlier than five days after peak. The lightcurves of 
these 18 supernovae have all been re -fit since P97, usin g the 
more recently available photometry ( frfamuy et al. 199*6 ) and 
our K corrections. 

Figures |l| and ||(a) show the Hubble diagram of effective rest- 
frame B magnitude corrected for the width-luminosity relation, 



effective 



= m R + A con -K B 



(3) 



as a function of redshift for the 42 Supernova Cosmol- 
ogy Project high-redshift supernovae, along with the 18 
Calan/Tololo low-redshift supernovae. (Here, K B r is the cross- 
filter K correction from observed R band to restframe B band.) 
Tables |l] and ^] give the corresponding IAU names, redshifts, 
magnitudes, corrected magnitudes, and their respective uncer- 
tainties. As inP97, the inner error bars in Figures |l] and ^repre- 
sent the photometric uncertainty, while the outer error bars add 
in quadrature 0.17 magnitudes of intrinsic dispersion of SN la 
magnitudes that remain after applying the width-luminosity 
correction. For these plots, the slope of the width-brightness 
relation was taken to be a = 0.6, the best-fit value of Fit C dis- 
cussed below. (Since both the low- and high-redshift supernova 
light-curve widths are clustered rather closely around s = 1, as 
shown in Figure |], the exact choice of a does not change the 
Hubble diagram significantly.) The theoretical curves for a uni- 
verse with no cosmological constant are shown as solid lines, 
for a range of mass density, £1m = 0, 1,2. The dashed lines rep- 
resent alternative flat cosmologies, for which the total mass- 
energy density Om+^a = 1 (where Oa = A/3H^). The range of 
models shown are for (Qm.^a) = (0,1), (0.5,0.5), (1,0), which 
is covered by the matching solid line, and (1.5, -0.5). 

3. FITS TO C M AND O a 

The combined low- and high-redshift supernova datasets of 
Figure [l] are fit to the Friedman-Robertson- Walker magnitude- 
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redshift relation, expressed as in P97: 

^effective = mR + a(s _ 1) _ Kbr _ Ar (4) 

= M B + 5\ogV L (z;n M ,n A ), 

where T>i = HqcIi is the "Hubble-constant-free" luminosity dis- 
tance and Mb = Mb -5logHo + 25 is the "Hubble-constant- 
free" B-band absolute magnitude at maximum of a SN la with 
width 5=1. (These quantities are, respectively, calculated 
from theory or fit from apparent magnitudes and redshifts, both 
without any need for Hq. The cosmological-parameter results 
are thus also completely independent of Hq.) The details of 
the fitting procedure as presented in P97 were followed, ex- 
cept that both the low- and high-redshift supernovae were fit 
simultaneously, so that Mb and a, the slope of the width- 
luminosity relation, could also be fit in addition to the cos- 
mological parameters J7m and For most of the analy- 

ses in this paper, Mb and a are statistical "nuisance" pa- 
rameters; we calculate 2-dimensional confidence regions and 
single-parameter uncertainties for the cosmological parame- 
ters by integrating over these parameters, i.e., V(Hm,^a) = 
J J V(QM,^A,M B ,a)dMB da. 

As in P97, the small correlations between the photometric 
uncertainties of the high-redshift supernovae, due to shared cal- 
ibration data, have been accounted for by fitting with a correla- 
tio n matrix of uncertainties. (The correlation matrix is available 
at frttp://www-supernova.rbl.gov .) The low-redshift supernova 
photometry is more likely to be uncorrelated in its calibration 
since these supernovae were not discovered in batches. How- 
ever, we take a 0.01 mag systematic uncertainty in the com- 
parison of the low-redshift B-band photometry and the high- 
redshift B-band photometry. The stretch-factor uncertainty is 
propagated with a fixed width-luminosity slope (taken from the 
low-redshift supernovae; cf. P97), and checked for consistency 
after the fit. 

We have compared the results of Bayesian and classical, "fre- 
quentist," fitting procedures. For the Bayesian fits, we have 
assumed a "prior" probability distribution that has zero proba- 
bility for f^M < 0, but otherwise uniform probability in the four 
parameters Hm, &a, a, and Mb- For the frequentist fits, we 
have followed the classical statistical procedures described by 



Feldman & Cousins (1998), to guarantee frequentist coverage 



of our confidence regions in the physically allowed part of pa- 
rameter space. Note that throughout the previous cosmology lit- 
erature, completely unconstrained fits have generally been used 
that can (and do) lead to confidence regions that include the part 
of parameter space with negative values for £7m- The differ- 
ences between the confidence regions that result from Bayesian 
and classical analyses are small. We present the Bayesian con- 
fidence regions in the figures, since they are somewhat more 
conservative, i.e. have larger confidence regions, in the vicinity 
of particular interest near A = 0. 

The residual dispersion in SN la peak magnitude after cor- 
recting for the width-luminosity relation is small, about 0.17 
magnitudes, before applying any color-correction. This was 
reported in Hamuy et al. (1996) for the low-redshift Calan- 
Tololo supernovae, and it is striking that the same residual is 
most consistent with the current 42 high-redshift supernovae 
(see Section 5). It is not clear from the current datasets, how- 
ever, whether this dispersion is best modeled as a normal dis- 
tribution (a Gaussian in flux space) or a log-normal distribution 
(a Gaussian in magnitude space). We have therefore performed 
the fits two ways: minimizing x 2 measured using either mag- 



nitude residuals or flux residuals. The results are generally in 
excellent agreement, but since the magnitude fits yield slightly 
larger confidence regions, we have again chosen this more con- 
servative alternative to report in this paper. 

We have analyzed the total set of 60 low- plus high-redshift 
supernovae in several ways, with the results of each fit pre- 
sented as a row of Table The most inclusive analyses are 
presented in the first two rows: Fit A is a fit to the entire dataset, 
while Fit B excludes two supernovae that are the most sig- 
nificant outliers from the average lightcurve width, 5=1, and 
two of the remaining supernovae that are the largest residu- 
als from Fit A. Figure Q shows that the remaining low- and 
high-redshift supernovae are well matched in their lightcurve 
width — the error-weighted means are (5)e a muy = 0.99 ± 0.01 
and (5)scp = 1.00 ±0.01 — making the results robust with re- 
spect to the width-luminosity-relation correction (see Section 
4.5). Our primary analysis, Fit C, further excludes two super- 
novae that are likely to be reddened, and is discussed in the 
following section. 

Fits A and B give very similar results. Removing the two 
large-residual supernovae from Fit A yields indistinguishable 
results, while Figure |^(a) shows that the 68% and 90% joint 
confidence regions for J1m and £!a still change very little af- 
ter also removing the two supernovae with outlier lightcurve 
widths. The best-fit mass-density in a flat universe for Fit A 
is, within a fraction of the uncertainty, the same value as for 



Fit b, n 



flat . 
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(see Table Bh. The main difference be- 



tween the fits is the goodness-of-fit: the larger x per degree of 
freedom for Fit A, xt = 1 -76, indicates that the outlier super- 
novae included in this fit are probably not part of a Gaussian 
distribution and thus will not be appropriately weighted in a x 2 
fit. The x 2 P er degree of freedom for Fit B, xt = 116, is over 
300 times more probable than that of fit A, and indicates that 
the remaining 56 supernovae are a reasonable fit to the model, 
with no large statistical errors remaining unaccounted for. 

Of the two large-residual supernovae excluded from the fits 
after Fit A, one is fainter than the best-fit prediction and one 
is brighter. The photometric color excess (see Section 4.1) for 
the fainter supernova, SN 19970, has an uncertainty that is too 
large to determine conclusively whether it is reddened. The 
brighter supernova, SN 1994H, is one of the first seven high- 
redshift supernovae originally analyzed in P97, and is one of 
the few supernovae without a spectrum to confirm its classifi- 
cation as a SN la. After re-analysis with additional calibration 
data and improved /f-corrections, it remains the brightest out- 
lier in the current sample, but it affects the final cosmological 
fits much less as one of 42 supernovae, rather than 1 of 5 super- 
novae in the primary P97 analysis. 



4. SYSTEMATIC UNCERTAINTIES AND CROSS-CHECKS 

With our large sample of 42 high-redshift SNe, it is not only 
possible to obtain good statistical uncertainties on the measured 
parameters, but also to quantify several possible sources of sys- 
tematic uncertainties. As discussed in P97, the primary ap- 
proach is to examine subsets of our data that will be affected 
to lesser extents by the systematic uncertainty being consid- 
ered. The high-redshift sample is now large enough that these 
subsets each contain enough supernovae to yield results of high 
statistical significance. 
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4. 1 . Extragalactic Extinction. 
4.1.1. Color-Excess Distributions 

Although we have accounted for extinction due to our 
Galaxy, it is still probable that some supernovae are dimmed by 
host galaxy dust or intergalactic dust. For a st andard dust ex- 
tinction law ( |Cardelli, Clayton, & Mathis 1989| ) the color, B-V, 
of a supernova will become redder as the amount of extinction, 
Ab, increases. We thus can look for any extinction differences 
between the low- and high-redshift supernovae by comparing 
their restframe colors. Since there is a small dependence of 
intrinsic color on the lightcurve width, supernova colors can 
only be compared for the same stretch factor; for a more conve- 
nient analysis, we subtract out the intrinsic colors, so that the re- 
maining color excesses can be compared simultaneously for all 
stretch factors. To determine the restframe color excess E(B-V) 
for each supernova, we fit the rest-frame B and V photometry to 
the B and V SN la lightcurve templates, with one of the fitting 
parameters representing the magnitude difference between the 
two bands at their respective peaks. Note that these lightcurve 
peaks are —2 days apart, so the resulting B lmx -V mdx color pa- 
rameter, which is frequently used to describe supernova colors, 
is not a color measurement on a particular day. The difference 
of this color parameter from the B m . dx -V max found for a sam- 
ple of l ow-redshift supernovae for the same lightcu rve stretch- 
factor ( fTripp 1998| ; [Kim et al. 1998) ; |Phillips 1998] ) does yield 
the restframe E(B-V) color excess for the fitted supernova. 

For the high-redshift supernovae at 0.3 < z < 0.7, the match- 
ing R- and /-band measurements take the place of the rest- 
frame B and V measurements and the fit B and V lightcurve 
templates are /^-corrected f rom the appropriate matching fil 
ters, e.g. R(t) = B(t) +K BR {t) (|Kim, Goobar, & Perlmutter 1996 



Nugent et al. 1998). For the three supernovae at z > 0.75, the 



observed R-I corresponds more closely to a restframe U-B 
color than to a B-V color, so E(B-V) is calculated from rest- 
frame E(U-B) using the extinction law of C ardelli, Clayton, 
& Mathis (1989). Similarly, for the two SNe la at z ~ 0.18, 
E(B-V) is calculated from restframe E(V-R). 

Figure ^ shows the color excess distributions for both the 
low- and high-redshift supernovae, after removing the color ex- 
cess due to our Galaxy. Six high-redshift supernovae are not 
shown on this E(B-V) plot, because six of the first seven high- 
redshift supernovae discovered were not observed in both R and 
/ bands. The color of one low-redshift supernova, SN 1992bc, 
is poorly determined by the V-band template fit and has also 
been excluded. Two supernovae in the high-redshift sample are 
> 3a red-and-faint outliers from the mean in the joint proba- 
bility distribution of E(B-V) color excess and magnitude resid- 
ual from Fit B. These two, SNe 1996cg and 1996cn (shown in 
light shading in Figure ||), are very likely reddened supernovae. 
To obtain a more robust fit of the cosmological parameters, in 
Fit C we remove these supernovae from the sample. As can 
be seen from the Fit-C 68% confidence region of Figure |5](a), 
these likely -reddened supernovae do not significantly affect any 
of our results. The main distribution of 38 high-redshift super- 
novae thus is barely affected by a few reddened events. We find 
identical results if we exclude the six supernovae without color 
measurements (Fit G in Table ||). We take Fit C to be our pri- 
mary analysis for this paper, and in Figure ^|, we show a more 
extensive range of confidence regions for this fit. 



4.1.2. Cross-checks on Extinction 

The color-excess distributions of the Fit C dataset (with the 
most significant measurements highlighted by dark shading in 
Figure show no significant difference between the low- and 
high-redshift means. The dashed curve drawn over the high- 
redshift distribution of Figure || shows the expected distribution 
if the low-redshift distribution had the measurement uncertain- 
ties of the high-redshift supernovae indicated by the dark shad- 
ing. This shows that the reddening distribution for the high- 
redshift SNe is consistent with the reddening distribution for the 
low-redshift SNe, within the measurement uncertainties. The 
error-weighted means of the low- and high-redshift distribu- 
tions are almost identical: (£'(-S-V))Hamuy = 0.033 ± 0.014 mag 
and (E{B-V)) SC p = 0.035 ±0.022 mag. We also find no sig- 
nificant correlation between the color excess and the statistical 
weight or redshift of the supernovae within these two redshift 
ranges. 

To test the effect of any remaining high-redshift reddening 
on the Fit C measurement of the cosmological parameters, we 
have constructed a Fit H-subset of the high-redshift supernovae 
that is intentional biased to be bluer than the low-redshift sam- 
ple. We exclude the error-weighted reddest 25% of the high- 
redshift supernovae; this excludes 9 high-redshift supernovae 
with the highest error-weighted E(B-V). We further exclude 
two supernovae that have large uncertainties in E(B-V) but are 
significantly faint in their residual from Fit C. This is a some- 
what conservative cut since it removes the faintest of the high- 
redshift supernovae, but it does ensure that the error-weighted 
E(B-V) mean of the remaining supernova subset is a good 
indicator of any reddening that could affect the cosmological 
parameters. The probability that the high-redshift subset of 
Fit H is redder in the mean than the low-redshift supernovae 
is less than 5%; This subset is thus very unlikely to be biased to 
fainter magnitudes by high-redshift reddening. Even with non- 
standard, "greyer" dust that does not cause as much reddening 
for the same amount of extinction, a conservative estimate of 
the probability that the high-redshift subset of Fit H is redder 
in the mean than the low-redshift supernovae is still less than 
-17%, for any high-redshift value of TZ B = A B /E(B-V) less 
than twice the low-redshift value. (These same confidence lev- 
els are obtained whether using Gaussian statistics, assuming a 
normal distribution of E(B-V) measurements, or using boot- 
strap resampling statistics, based on the observed distribution.) 
The confidence regions of Figure |]^c) and the il^ results in 
Table || show that the cosmological parameters found for Fit H 
differ by less than half of a standard deviation from those for 
Fit C. We take the difference of these fits, 0.03 in ft£f (which 
corresponds to less than 0.025 in magnitudes) as a — la upper 
bound on the systematic uncertainty due to extinction by dust 
that reddens. 

Note that the modes of both distributions appear to be at zero 
reddening, and similarly the medians of the distributions are 
quite close to zero reddening: (£(fi-V))]^™ = 0.01 mag and 

(E(B-V))f£f' dn = 0.00 mag. This should be taken as sugges- 
tive rather than conclusive since the zeropoint of the relation- 
ship between true color and stretch is not tightly constrained 
by the current low-redshift SN la dataset. This apparent strong 
clustering of SNe la about zero reddening has been noted in 
the past for low-redshift supernova samples. Proposed expla- 
nations have been given based on the relative spatial distribu- 
tions of the SNe la and the dust: Modeling by H atano, Branch, 
& Deaton (1997) of the expected extinction of SN la disk and 
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bulge populations viewed at random orientations shows an ex- 
tinction distribution with a strong spiked peak near zero ex- 
tinction along with a broad, lower-probability wing to higher 
extinction. This wing will be further suppressed by the obser- 
vational selection against more reddened SNe, since they are 
dimmer. (For a flux-limited survey this suppression factor is 

iQ-atWBEQWyals-T)] _ 1() -1.6£(iH0 where flR ifj the slope of the 

supernova number counts.) We also note that the high-redshift 
supernovae for which we have accurate measurements of appar- 
ent separation between SN and host position (generally, those 
with Hubble Space Telescope imaging) appear to be relatively 
far from the host center, despite our high sear ch sensitivity to 



supernovae in front of the host galax y core (see ^ ain et al 



for se arch efficiency studies; also cf. Wang, Htiflich, & Wheele 



1996 



1997). If generally true for the entire sample, this would be 



consistent with little extinction. 

Our results, however, do not depend on the low- and high- 
redshift color-excess distributions being consistent with zero 
reddening. It is only important that the reddening distributions 
for the low-redshift and high-redshift datasets are statistically 
the same, and that there is no correlation between reddening 
and statistical weight in the fit of the cosmological parame- 
ters. With both of these conditions satisfied, we find that our 
measurement of the cosmological parameters is unaffected (to 
within the statistical error) by any small remaining extinction 
among the supernovae in the two datasets. 

4.1.3. Analysis with Reddening Correction of 
Individual Supernovae 

We have also performed fits using restframe fi-band mag- 
nitudes individually corrected for host galaxy extinction using 
Ab = IZbE(B-V) (implicitly assuming that the extragalactic ex- 
tinction is all at the redshift of the host galaxy). As a direct 
comparison between the treatment of host galaxy extinction de- 
scribed above an d an alternative Bayesian method (R iess, Press, 
& Kirshner 1996), we applied it to the 53 SNe la with color 
measurements in our Fit C dataset. We find that our cosmo- 
logical parameter results are robust with respect to this change, 
although this method can introduce a bias into the extinction 
corrections, and hence the cosmological parameters. In brief, 
in this method the Gaussian extinction probability distribution 
implied by the measured color-excess and its error is multiplied 
by an assumed a priori probability distribution (the Bayesian 
prior) for the intrinsic distribution of host extinctions. The 
most probable value of the resulting renormalized probability 
distribution is taken as the extinction, and following Riess (pri- 
vate communication) the second-moment is taken as the uncer- 
tainty. For this analysis, we choose a conservative prior (as 
given in Riess, Press, & Kirshner 1996) that does not assume 



that the supernovae are unextinguished, but rather is somewhat 
broader than the true extinction distribution where the majority 
of the previously observed supernovae apparently suffer very 
little reddening. (If one alternatively assumes that the current 
data's extinction distribution is quite as narrow as that of previ- 
ously observed supernovae, one can choose a less conservative 
but mo re realistic narrow prior probability distribution, such as 
that of Hatano, Branch, & Deaton (1997). This turns out to 



be quite similar to our pre vious analysis in Section 4. 1 . 1 , since 
a distribution like that of Hatano, Branch, & Deaton has zero 



extinction for most supernovae.) 

This Bayesian method with a conservative prior will only 
brighten supernovae, never make them fainter, since it only af- 
fects the supernovae with redder measurements than the zero- 



extinction E(B—V) value, leaving unchanged those measured to 
be bluer than this. The resulting slight difference between the 
assumed and true reddening distributions would make no differ- 
ence in the cosmology measurements if its size were the same 
at low and high redshifts. However, since the uncertainties, 
CT E(zw)> m me high-redshift dataset E(B-V) measurements are 
larger on average than those of the low-redshift dataset, cr^^, 
this method can over-correct the high-redshift supernovae on 
average relative to the low-redshift supernovae. Fortunately, as 
shown in Appendix A, even an extreme case with a true dis- 
tribution all at zero extinction and a conservative prior would 
introduce a bias in extinction Ab only of order 0.1 magnitudes 
at worst for our current low- and high-redshift measurement un- 
certainties. The results of Fit E are shown in Table || and as the 
dashed contour in Figure |]^d), where it can be seen that com- 
pared to Fit C this approach moves the best fit value much less 
than this, and in the direction expected for this effect (indicated 
by the arrows in Figure |5]d). The fact that 57^" changes so little 
from Case C, even with the possible bias, gives further confi- 
dence in the cosmological results. 

We can eliminate any such small bias of this method by as- 
suming no Bayesian prior on the host-galaxy extinction, allow- 
ing extinction corrections to be negative in the case of super- 
novae measured to be bluer than the zero-extinction E(B-V) 
value. As expected, we recover the unbiased results within er- 
ror, but with larger uncertainties since the Bayesian prior also 
narrows t he error bars in the method of R iess, Press, & Kirsh- 
ner (1996). However, there remains a potential source of bias 
when correcting for reddening: the effective ratio of total to se- 
lective extinction, TZb, could vary, for several reasons. First, 
the extinction could be due to host galaxy dust at the super- 
nova's redshift or intergalactic dust at lower redshifts, where it 
will redden the supernova less since it is acting on a redshifted 
spectrum. Second, TZb may be sensitive to dust density, as indi- 
cated by variations in the dust extinction laws betw een various 
sight-lines in the Galaxy (|Clayton & Cardelli 1988|; Gordon & 
Clayton 1998). Changes in metallicity might be expected to 
be a third possible cause of TZb evolution, since metallicity is 
one dust-relate d quantity known to evolve with redshift (P et- 
tini et al. 1997), but fortunately it appears not to significantly 
alter TZb as evidenced by the similarity of the optical portions 
of the extinction curves of the Gala xy, the LMC, and the SMC 
( |Pei 1992| ; |Gordon & Clayton 1998[ ). Three-filter photometry of 
high-redshift supernovae currently in progress with the Hubble 
Space Telescope will help test for such differences in TZb- 

To avoid these sources of bias, we consider it important to 
use and compare both analysis approaches: the rejection of red- 
dened supernovae and the correction of reddened supernovae. 
We do find consistency in the results calculated both ways. The 
advantages of the analyses with reddening corrections applied 
to individual supernovae (with or without a Bayesian prior on 
host-galaxy extinction) are outweighed by the disadvantages for 
our sample of high-redshift supernovae; although, in principle, 
by applying reddening corrections the intrinsic magnitude dis- 
persion of SNe la can be reduced from an observed dispersion 
of 0.17 magnitudes to approximately 0.12 magnitudes, in prac- 
tice the net improvement for our sample is not significant since 
uncertainties in the color measurements often dominate. We 
have therefore chosen for our primary analysis to follow the first 
procedure discussed above, removing the likely-reddened su- 
pernovae (Fit C) and then comparing color-excess means. The 
systematic difference for Fit H, which rejects the reddest and 
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the faintest high-redshift supernovae, is already quite small, and 
we avoid introducing additional actual and possible biases. Of 
course, neither approach avoids biases if TZb at high redshift 
is so large [> 21Zb(z = 0)] that dust does not redden the su- 
pernovae enough to be distinguished and this dust makes more 
than a few supernovae faint. 

4.2. Malmquist Bias and other Luminosity Biases. 

In the fit of the cosmological parameters to the magnitude- 
redshift relation, the low-redshift supernova magnitudes pri- 
marily determine and the width-luminosity slope a, and 
then the comparison with the high-redshift supernova magni- 
tudes primarily determines £1m and Qjy, Both low- and high- 
redshift supernova samples can be biased towards selecting the 
brighter tail of any distribution in supernova detection magni- 
tude for supernovae found near th e detection thres h old of the 
search (classical Malmquist bias; ^vlalmquist 1924, 1936). A 
width-luminosity relation fit to such a biased population would 
have a slope that is slightly too shallow and a zeropoint slightly 
too bright. A second bias is also acting on the supernova sam- 
ples, selecting against supernovae on the narrow-lightcurve side 
of the width-luminosity relation since such supernovae are de- 
tectable for a shorter period of time. Since this bias removes the 
narrowest/faintest supernova lightcurves preferentially, it culls 
out the part of the width-brightness distribution most subject 
to Malmquist bias, and moves the resulting best-fit slope and 
zeropoint closer to their correct values. 

If the Malmquist bias is the same in both datasets, then it is 
completely absorbed by A4b and a and does not affect the cos- 
mological parameters. Thus, our principal concern is that there 
could be a difference in the amount of bias between the low- 
redshift and high-redshift samples. Note that effects peculiar to 
photographic SNe searches, such as saturation in galaxy cores, 
which might in principle select slightly different SNe la sub- 
populations should not be important in determining luminosity 
bias because lightcurve stretch compensates for any such differ- 
ences. Moreover, Figure Q shows that the high-redshift SNe la 
we have discovered have a stretch distribution entirely consis- 
tent with those discovered in the Calan/Tololo search. 

To estimate the Malmquist bias of the high-redshift- 
supernova sample, we first determined the completeness of our 
high-redshift searches as a function of magnitude, through an 
exte nsive series of tes ts inserting artificial SNe into our images 
(see ^ain et al. 1996 ). We find that roughly 30% of our high- 
redshift supernovae were detected within twice the SN la intrin- 
sic luminosity dispersion of the 50% completeness limit, where 
the above biases might be important. This is consistent with 
a simple model where the supernova number counts follow a 
power-law slope of .4 mag" 1 , similar t o that seen for compara- 
bly distant galaxies ( |Smail et al. 1995 ). For a flux-limited sur- 
vey of standard candles having the lightcurve- width-corrected 
luminosity dispersion for SN la of ^0. 17 mag and this number- 
count power-law slope, we can calculate that the classical 
Malmqui st bias should be 0.03 mag (see, e.g., Mihalas & Bin- 
ney 1981, for a derivation of the classical Malmquist bias). 
(Note that this estimate is much smaller than the Malmquist 
bias affecting other cosmological distance indicators, due to the 
much smaller intrinsic luminosity dispersion of SNe la.) These 
high-redshift supernovae, however, are typically detected be- 
fore maximum, and their detection magnitudes and peak mag- 
nitudes have a correlation coefficient of only 0.35, so the ef- 
fects of classical Malmquist bi as should be diluted. Applying 
the formalism of Willick (1994) we estimate that the decorrela- 



tion between detection magnitude and peak magnitude reduces 
the classical Malmquist bias in the high-redshift sample to only 

0. 01 mag. The redshift and stretch distributions of the high- 
redshift supernovae that are near the 50%-completeness limit 
track those of the overall high-redshift sample, again suggest- 
ing that Malmquist biases are small for our dataset. 

We cannot make an exactly parallel estimate of Malmquist 
bias for the low-redshift-supernova sample, because we do not 
have information for the Calan/Tololo dataset concerning the 
number of supernovae found near the detection limit. However, 
the amount of classical Malmquist bias should be similar for 
the Calan/Tololo SNe since the amount of bias is dominated by 
the intrinsic luminosity dispersion of SNe la, which we find to 
be the same for the low-redshift and high-redshift samples (see 
Section 5). Figure ^ shows that the stretch distributions for the 
high-redshift and low-redshift samples are very similar, so that 
the compensating effects of stretch-bias should also be similar 
in the two datasets. The major source of difference in the bias 
is expected to be due to the close correlation between the de- 
tection magnitude and the peak magnitude for the low-redshift 
supernova search, since this search tended not to find the super- 
novae as early before peak as the high-redshift search. In ad- 
dition, t he number-counts at low-redshift should be somewhat 
steeper ( [Vladdox et al. 1990 ). We thus expect the Calan/Tololo 
SNe to have a bias closer to that obtained by direct application 
of the the classical Malmquist bias formula, 0.04 mag. One 
might also expect "inhomogeneous Malmquist bias" to be more 
important for the low-redshift supernovae, since in smaller vol- 
umes of space inhomogeneities in the host galaxy distribution 
might by chance put more supernovae near the detection limit 
than would be expected for a homogeneous distribution. How- 
ever, after averaging over all the Calan/Tololo supernova-search 
fields the total low-redshift volume searched is large enough 
that we expect galaxy count fluctuations of only ~4%, so the 
classical Malmquist bias is still a good approximation. 

We believe that both these low- and high-redshift biases may 
be smaller, and even closer to each other, due to the mitigating 
effect of the bias against detection of low-stretch supernovae, 
discussed above. However, to be conservative, we take the clas- 
sical Malmquist bias of 0.04 mag for the low-redshift dataset, 
and the least biased value of 0.01 mag for the high-redshift 
dataset, and consider systematic uncertainty from this source to 
be the difference, 0.03 mag, in the direction of low-redshift su- 
pernovae more biased than high-redshift. In the other direction, 

1. e. for high-redshift supernovae more biased than low-redshift, 
we consider the extreme case of a fortuitously unbiased low- 
redshift sample, and take the systematic uncertainty bound to 
be the 0.01 mag bias of the high-redshift sample. (In this direc- 
tion any systematic error is less relevant to the question of the 
existence of a cosmological constant.) 

4.3. Gravitational Lensing. 

As discussed in P97, the clumping of mass in the universe 
could leave the line-of-sight to most of the supernovae under- 
dense, while occasional supernovae may be seen through over- 
dense regions. The latter supernovae could be significantly 
brightened by gravitational lensing, while the former super- 
novae would appear somewhat fainter. With enough super- 
novae, this effect will average out (for inclusive fits, such as 
Fit A, which include outliers), but the most over-dense lines 
of sight may be so rare that a set of 42 supernovae may only 
sample a slightly biased (fainter) set. The probability distribu- 
tion of these amplifications and deamplifications has previously 
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been studied both analytically and by Monte Carlo simulations. 
Given the acceptance window of our supernova search, we can 
integrate the probability distributions from these studies to es- 
timate the bias due to amplified or deamplified supernovae that 
may be rejected as outliers. This average (de)amplificationbias 
is less than 1% at the redshifts of our supernovae for simula- 
t ions based on isoth ermal spheres the size of typical galaxies 
( Holz & Wald 1998), N-body simulations using reali stic mass 
power spectra ( Wambsganss, Cen, & Ostriker 1998 ), and the 
analytic models of jFrieman (1996 ). 

It is also possible that the small-scale clumping of matter 
is more extreme, e.g., if significant amounts of mass were in 
the form of compact objects such as MACHOs. This could 
lead to many supernova sightlines that are not just under-dense, 
but nearly empty. Once again, only the very rare line of sight 
would have a compact object in it, amplifying the supernova 
signal. To first approximation, with 42 supernovae we would 
see only the nearly empty beams, and thus only deamplifi- 
cations. The appropriate luminosity-distance formula in this 
case is not the Friedmann-Robertson- Walker (FRW) formula 
but rather the "partia lly filled beam" f ormula with a mass fill- 
ing factor, rj w (see Kantowski 1998, and references therein). 



We present the results of the fit of our data (Fit K) with this 
luminosity-distance formula (as ca lculated using the code of 
Kayser, Helbig, & Schramm 1996) in Figure pT A more re- 



alistic limit on this point-like mass density can be estimated, 
because we would expect such point-like masses to collect into 
the gravi tational potential wells already mark ed by galaxies and 
clusters, ^ukugita, Hogan, & Peebles (1997 ) estimate an upper 
limit of J1m < 0.25 on the mass which is clustered like galaxies. 
In Figure pi we also show the confidence region from Fit L, as- 
suming that only the mass density contribution up to Qm = 0.25 
is point-like, with filling factor r] = 0, and that r\ rises to 0.75 
at £!m = 1. We see that at low mass density, the Friedman- 
Robertson- Walker fit is already very close to the nearly empty- 
beam (rj w 0) scenario, so the results are quite similar. At high 
mass density, the results diverge, although only minimally for 



Fit L; the best fit in a flat universe is fijjf : 



o 4 +().io 
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4.4. Supernova Evolution and Progenitor Environment 
Evolution 

The spectrum of a SN la on any given point in its lightcurve 
reflects the complex physical state of the supernova on that 
day: the distribution, abundances, excitations, and velocities 
of the elements that the photons encounter as they leave the 
expanding photosphere all imprint on the spectra. So far, the 
high-redshift supernovae that have been studied have lightcurve 
shapes just like those of low-redshift supernovae (see G old- 
haber et al. 1998), and their spectra show the same features on 
the same day of the lightcurve as their low-redshift counter- 
parts having comparable lightcurve width. This is true all the 
way out to the z = 0.83 limit of the current sample ( Perlmutter 
et al. 1998b). We take this as a strong indication that the phys- 
ical parameters of the supernova explosions are not evolving 
significantly over this time span. 

Theoretically, evolutionary effects might be caused by 
changes in progenitor populations or environments. For ex- 
ample, lower metallicity and more massive SN la-progenitor 
binary systems should be found in younger stellar populations. 
For the redshifts that we are considering, z < 0.85, the change 
in average progenitor masses may be small (R uiz-Lapuente, 
Canal, & Burkert 1997; Ruiz-Lapuente 1998). However, such 
progenitor mass differences or differences in typical progenitor 



metallicity are expected to lead to differences in the final C/O 
ratio in the exploding white dwarf, and hence affect the energet- 
ics of the explosion. The primary concern here would be if this 
changed the zero-point of the width-luminosity relation. We 
can look for such changes by comparing lightcurve rise times 
between low and high-redshift supernova samples, since this is 
a sensitive indicator of explosion energetics. Preliminary in- 
dications suggest that no significant rise-time change is seen, 
with an upper limit of < 1 day for our sample (see forthcoming 
high-redshift studies of Goldhaber e t al. 1998 and Nugent et al . 
1998, and low-redshift b oun ds from Vacca & Leibundgut 1996, 
Leibundgut et al.| |l996b[ and Marvin & Perlmutter ||1 9 8 9| ). This 
tight a constraint on rise-time change would theoretically limit 
the z ero-point change to less than ^0. 1 mag (see N ugent et al. 
1995; Hoflich, Wheeler, & Thielemann 199ty . 

A change in typical C/O ratio can also affect the ignition 
density of the explosion and the propagation characteristics of 
the burning front. Such changes would be expected to appear 
as diff erences in lightcurve timesc ales before and after max- 
imum (Hoflich & Khokhlov 1996). Preliminary indications 



of consistency between such low- and high-redshift lightcurve 
timescales suggest that this is probably not a major effect for 
our supernova samples (Goldhaber et al., 1998). 

Changes in typical progenitor metallicity should also directly 
cause some differences in SN la spectral features (H oflich, 
Wheeler, & Thielemann 1998). Spectral differences big enough 
to affect the B and V-band lightcurves (see, for example, the ex- 
treme mixing models presented in Figure 9 of H oflich, Wheeler, 
& Thielemann 1998) should be clearly visible for the best 
signal-to-noise spectra we h ave obtained for our di stant super- 
novae, yet t hey ar e not seen ( [Filippenko et al. 1998 ; Hook, Nu- 
gent, et al., 1998). The consistency of slopes in the lightcurve 
width-luminosity relation for the low- and high-redshift super- 
novae can also constr ain the possibility of a strong metallicity 
effect of the type that frloflich, Wheeler, & Thielemann (1998 ) 
describes. 

An additional concern might be that even small changes in 
spectral features with metallicity could in turn affect the calcu- 
lations of K corrections and reddening corrections. This effect, 
too, is very small, less than 0.01 magnitudes, for photomet- 
ric observations of SNe la conducted in the restframe B or V 
bands (see Figures 8 and 10 of Hoflich, Wheeler, & Thiele- 
mann 1998), as is the case for almost all of our supernovae. 
(Only two of our supernovae have primary observations that 
are sensitive to the restframe U band, where the magnitude can 
change by ^0.05 magnitudes, and these are the two supernovae 
with the lowest weights in our fits, as shown by the error bars 
of Figures ^| In general the /-band observations, which are 
mostly sensitive to the restframe B band, provide the primary 
lightcurve at redshifts above 0.7.) 

The above analyses constrain only the effect of progenitor- 
environment evolution on SN la intrinsic luminosity; however, 
the extinction of the supernova light could also be affected, 
if the amount or character of the dust evolves, e.g. with host 
galaxy age. In Section 4.1, we limited the size of this extinc- 
tion evolution for dust that reddens, but evolution of "grey" 
dust grains larger than ^0. 1 /im, which would cause more color- 
neutral optical extinction, can evade these color measurements. 
The following two analysis approaches can constrain both evo- 
lution effects, intrinsic SN la luminosity evolution and extinc- 
tion evolution. They take advantage of the fact that galaxy 
properties such as formation age, star-formation history, and 
metallicity are not monotonic functions of redshift, so even the 
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low-redshift SNe la are found in galaxies with a wide range of 
ages and metallicities. It is a shift in the distribution of relevant 
host-galaxy properties occurring between z ~ and z ~ 0.5 that 
could cause any evolutionary effects. 

Width-Luminosity Relation Across Low-Redshift Environ- 
ments. To the extent that low-redshift SNe la arise from pro- 
genitors with a range of metallicities and ages, the lightcurve 
width-luminosity relation dis covered for these SN e can already 
account for these effects (cf. |Hamuy et al. 1995 , 1996). When 
corrected for the width-luminosity relation, the peak magni- 
tudes of low-redshift SNe la exhibit a very narrow magnitude 
dispersion about the Hubble line, with no evidence of a signif- 
icant progenitor-environment difference in the residuals from 
this fit. It therefore does not matter if the population of progen- 
itors evolves such that the measured lightcurve widths change, 
since the width-luminosity relation apparently is able to correct 
for these changes. It will be important to continue to study fur- 
ther nearby SNe la to test this conclusion with as wide a range 
of host-galaxy ages and metallicities as possible. 

Matching Low- and High-Redshift Environments. Galaxies 
with different morphological classifications result from differ- 
ent evolutionary histories. To the extent that galaxies with sim- 
ilar classifications have similar histories, we can also check 
for evolutionary effects by using supernovae in our cosmol- 
ogy measurements with matching host galaxy classifications. 
If the same cosmological results are found for each measure- 
ment based on a subset of low- and high-redshift supernovae 
sharing a given host-galaxy classification, we can rule out many 
evolutionary scenarios. In the simplest such test, we compare 
the cosmological parameters measured from low- and high- 
redshift elliptical host galaxies with those measured from low- 
and high-redshift spiral host galaxies. Without high-resolution 
host-galaxy images for most of our high-redshift sample, we 
currently can only approximate this test for the smaller num- 
ber of supernovae for which the host-galaxy spectrum gives a 
strong indication of galaxy classification. The resulting sets of 
9 elliptical-host and 8 spiral-host high-redshift supernovae are 
matched to the 4 elliptical-host and 10 spiral-host low-redshift 
supernovae (based on the morphological classifications listed in 



Hamuy et al. 1996| , and excluding two with SBO hosts). We find 



no significant change in the best-fit cosmology for the elliptical 
host-galaxy subset (with both the low- and high-redshift subsets 
about one sigma brighter than the mean of the full sets), and a 
small ( < 1 er) shift lower in fijjf for the spiral host-galaxy subset. 
Although the consistency of these subset results is encouraging, 
the uncertainties are still large enough (approximately twice the 
Fit C uncertainties) that this test will need to await the host- 
galaxy classification of the full set of high-redshift supernovae 
and a larger low-redshift supernova sample. 

4.5. Further Cross-Checks 

We have checked several other possible effects that might 
bias our results, by fitting different supernova subsets and using 
alternative analyses: 

Sensitivity to Width-Luminosity Correction. Although the 
lightcurve width correction provides some insurance against su- 
pernova evolution biasing our results, Figure Q shows that al- 
most all of the Fit C supernovae at both low- and high-redshift 
are clustered tightly around the most-probable value of s = 1, 
the standard width for a B-band Leibundgut SN la template 
lightcurve. Our results are therefore rather robust with respect 



to the actual choice of width-luminosity relation. We have 
tested this sensitivity by re-fitting the supernovae of Fit C, but 
with no width-luminosity correction. The results (Fit D), as 
shown in Figure ||(b), and listed in Table ||, are in extremely 
close agreement with those of the lightcurve-width-corrected 
Fit C. The statistical uncertainties are also quite close; the 
lightcurve- width correction does not significantly improve the 
statistical dispersion for the magnitude residuals, because of the 
uncertainty in s, the measured lightcurve width. It is clear that 
the best-fit cosmology does not depend strongly on the extra 
degree of freedom allowed by including the width-luminosity 
relation in the fit. 

Sensitivity to Non-SN la Contamination. We have tested 
for the possibility of contamination by non-SN la events mas- 
querading as SNe la in our sample, by performing a fit af- 
ter excluding any supernovae with less certain SN la spectro- 
scopic and photometric identification. This selection removes 
the large statistical outliers from the sample. In part, this may 
be because the host-galaxy contamination that can make it dif- 
ficult to identify the supernova spectrum can also increase the 
odds of extinction or other systematic uncertainties in photom- 
etry. For this more "pure" sample of 43 supernovae, we find 
Sljjf = 0.33±g go, just over half of a standard deviation from 
Fit C. 

Sensitivity to Galactic Extinction Model. Finally, we have 
tested the effect of the choi ce of Galactic extinction model, 
with a fit using the model of burstein & Heiles (1982| ), rather 
than Schlegel, Finkbeiner, & Davis (1998 ). We find no sig- 
nificant difference in the best-fit cosmological parameters, al- 
though we note that th e extinction near the Galactic p ole is 
somewhat larger in the Schlegel, Finkbeiner, & Davis model 
and this leads to a ^0.03 magnitude larger average offset be- 
tween the low-redshift supernova B-band observations and the 
high-redshift supernovae B-band observations. 

5. RESULTS AND ERROR BUDGET 

From Table || and Figure ||(a), it is clear that the results of 
Fits A, B, and C are quite close to each other, so we can con- 
clude that our measurement is robust with respect to the choice 
of these supernova subsets. The inclusive Fits A and B are the 
fits with the least subjective selection of the data. They already 
indicate the main cosmological results from this dataset. How- 
ever, to make our results robust with respect to host-galaxy 
reddening, we use Fit C as our primary fit in this paper. For 
Fit C, we find fijjf = 0.28^ in a flat universe. Cosmolo- 
gies with Qa = are a poor fit to the data, at the 99.8% con- 
fidence level. The contours of Figure [7] more fully charac- 
terize the best-fit confidence regions. (The table of this two- 
dimensional proba bility distribution is available at h ttp://www- 
supernova.lbl.gov/.) 

The residual plots of Figure ^(b and c) indicate that the best- 
fit f^ 1 in a flat universe is consistent across the redshift range 
of the high-redshift supernovae. Figure ||(c) shows the residu- 
als normalized by uncertainties; their scatter can be seen to be 
typical of a normal-distributed variable, with the exception of 
the two outlier supernovae that are removed from all fits after 
Fit A, as discussed above. Figure || compares the magnitude- 
residual distributions (the projections of Figure ||b) to the Gaus- 
sian distributions expected given the measurement uncertain- 
ties and an intrinsic dispersion of 0.17 mag. Both the low- 
and high-redshift distributions are consistent with the expected 
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distributions; the formal calculation of the SN la intrinsic- 
dispersion component of the observed magnitude dispersion 

) yields extrinsic =0.154 ±0.04 



observed 



measurement 



(it 

v intrinsi 

for the low-redshift distribution and c intrinsic 

= 0.157 ±0.025 for 

the high-redshift distribution. The % 2 per degree of freedom for 
this fit, xt = 1-12, also indicates that the fit model is a reason- 
able description of the data. The narrow intrinsic dispersion — 
which does not increase at high redshift — provides additional 
evidence against an increase in extinction with redshift. Even 
if there is grey dust that dims the supernovae without redden- 
ing them, the dispersion would increase, unless the dust is dis- 
tributed very uniformly. 

A flat, £1a = cosmology is a quite poor fit to the data. 
The (17m, ^a) = (1)0) line on Figure ||(b) shows that 38 out of 
42 high-redshift supernovae are fainter than predicted for this 
model. These supernovae would have to be over 0.4 magni- 
tudes brighter than measured (or the low-redshift supernovae 
0.4 magnitudes fainter) for this model to fit the data. 

The (£1m,^a) = (0,0) upper solid line on Figure ^](a) shows 
that the data are still not a good fit to an "empty universe," with 
zero mass density and cosmological constant. The high-redshift 
supernovae are as a group fainter than predicted for this cos- 
mology; in this case, these supernovae would have to be almost 
0.15 magnitudes brighter for this empty cosmology to fit the 
data, and the discrepancy is even larger for 17m > 0. This is 
reflected in the hig h probability (99.8%) of Q\ > 0. 

As discussed in 5oob_ar & Perlmutter (1995 ), the slope of 
the contours in Figure [7] is a function of the supernova red- 
shift distribution; since most of the supernovae reported here 
are near z ~ 0.5, the confidence region is approximately fit by 
0.8 £1m — 0.617a w -0.2 ±0.1. (The orthogonal linear combi- 
nation, which is poorly constrained, is fit by 0.6 17m + 0.8 S7a ~ 
1.5 ±0.7.) In P97, we emphasized that the well-constrained 
linear combination is not parallel to any contour of constant 
current-deceleration-parameter, qo = 17m/2— 17 a! the accelerat- 
ing/decelerating universe line of Figure ^| shows one such con- 
tour at qo = 0. Note that with almost all of the confidence region 
above this line, only currently accelerating universes fit the data 
well. As more of our highest redshift supernovae are analyzed, 
the long dimension of the confidence region will shorten. 

Error Budget 

Most of the sources of statistical error contribute a statistical 
uncertainty to each supernova individually, and are included in 
the uncertainties listed in Tables [l] and with small correla- 
tions between these u ncertainties given in the correla ted-error 
matrices (available at http://www-supernova.lbl.gov). These 



supernova-specific statistical uncertainties include the measure- 
ment errors on SN peak magnitude, lightcurve stretch factor, 
and absolute photometric calibration. The two sources of sta- 
tistical error that are common to all the supernovae are the in- 
trinsic dispersion of SN la luminosities after correcting for the 
width-luminosity relation, taken as 0.17 mag, and the redshift 
uncertainty due to peculiar velocities, which are taken as 300 
km s" 1 . Note that the statistical error in Mb and a are derived 
quantities from our four-parameter fits. By integrating the four- 
dimensional probability distributions over these two variables, 
their uncertainties are included in the final statistical errors. 

All uncertainties that are not included in the statistical error 
budget are treated as systematic errors for the purposes of this 
paper. In Sections 2 and 4, we have identified and bounded 
four potentially significant sources of systematic uncertainty: 
(1) the extinction uncertainty for dust that reddens, bounded 



at <0.025 magnitudes, the maximal effect of the nine reddest 
and two faintest of the high-redshift supernovae; (2) the differ- 
ence between the Malmquist bias of the low- and high-redshift 
supernovae, bounded at <0.03 magnitudes for low-redshift su- 
pernovae biased intrinsically brighter than high-redshift super- 
novae, and <0.01 magnitudes for high-redshift supernovae bi- 
ased brighter than low-redshift supernovae; (3) the cross-filter 
^-correction uncertainty of <0.02 magnitudes; and (4) the 
<0.01 magnitudes uncertainty in K corrections and reddening 
corrections due to the effect of progenitor metallicity evolution 
on the rest-frame B-band spectral features. We take the total 
identified systematic uncertainty to be the quadrature sum of 
the sources: +0.04 magnitudes in the direction of spuriously 
fainter high-redshift or brighter low-redshift supernovae, and 
-0.03 magnitudes in the opposite direction. 

Note that we treat the possibility of gravitational lensing due 
to small-scale clumping of mass as a separate analysis case, 
rather than as a contributing systematic error in our primary 
analysis; the total systematic uncertainty applies to this analy- 
sis as well. There are also several more hypothetical sources of 
systematic error discussed in Section 4, which are not included 
in our calculation of identified systematics. These include grey 
dust [with TZ B (z = 0.5) > 21Z B (z = 0)] and any SN la evolution- 
ary effects that might change the zero point of the lightcurve 
width-luminosity relation. We have presented bounds and tests 
for these effects which give preliminary indications that they 
are not large sources of uncertainty, but at this time they remain 
difficult to quantify, at least partly because the proposed phys- 
ical processes and entities that might cause the effects are not 
completely defined. 

To characterize the effect of the identified systematic uncer- 
tainties, we have refit the supernovae of Fit C for the hypothet- 
ical case (Fit J) in which each of the high-redshift supernovae 
were discovered to be 0.04 magnitudes brighter than measured, 
or, equivalently, the low-redshift supernovae were discovered 
to be 0.04 magnitudes fainter than measured. Figure ||(e) and 
Table [| show the results of this fit. The best-fit flat-universe 
Sl{Jf varies from that of Fit C by 0.05, less than the statistical 
error bar. The probability of 17a > is still over 99%. When 
we fit with the smaller systematic error in the opposite direc- 
tion (i.e., high-redshift supernovae discovered to be 0.03 mag- 
nitudes fainter than measured), we find (Fit I) only a 0.04 shift 
in 17£f from Fit C. 

The measurement error of the cosmological parameters has 
contributions from both the low- and high-redshift supernova 
datasets. To identify the approximate relative importance of 
these two contributory sources, we reanalyzed the Fit C dataset, 
first fitting Mb and a to the low-redshift dataset (this is rela- 
tively insensitive to cosmological model), and then fitting 17m 
and Ha to the high-redshift dataset. (This is only an approxi- 
mation, since it neglects the small influence of the low-redshift 
supernovae on S7m and S7a, and of the high-redshift supernovae 
on Mb and a, in the standard four-parameter fit.) Figure ||(f) 
shows this SIm-^a fit as a solid contour (labeled Fit M), with 
the 1-sigma uncertainties on Mb and a included with the sys- 
tematic uncertainties in the dashed-line confidence contours. 



This approach parallels the analyses of Perlmutter et al. ( 1 997e , 



1998c, 1998b), and thus also provides a direct comparison with 



the earlier results. We find that the more important contribution 
to the uncertainty is currently due to the low-redshift supernova 
sample. If three times as many well-observed low-redshift su- 
pernovae were discovered and included in the analysis, then the 
statistical uncertainty from the low-redshift dataset would be 
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smaller than the other sources of uncertainty. 

We summarize the relative statistical and systematic uncer- 
tainty contributions in Table f|. 

6. CONCLUSIONS AND DISCUSSION 

The confidence regions of Figure [7] and the residual plot of 
Figure ||(b) lead to several striking implications. First, the data 
are strongly inconsistent with the A = 0, flat universe model 
(indicated with a circle) that has been the theoretically favored 
cosmology. If the simplest inflationary theories are correct and 
the universe is spatially flat, then the supernova data imply that 
there is a significant, positive cosmological constant. Thus, the 
universe may be flat, or there may be little or no cosmological 
constant, but the data are not consistent with both possibilities 
simultaneously. This is the most unambiguous result of the cur- 
rent dataset. 

Second, this dataset directly addresses the age of the uni- 
verse relative to the Hubble time, Hq 1 , Figure ^ shows that 
the Hm—^Ia confidence regions are almost parallel to contours 
of constant age. For any value of the Hubble constant less 
than Ho = 70 km s" 1 Mpc -1 , the implied age of the universe 
is greater than 13 Gyr, allowing enough time for the old est stars 
in globula r clusters to evolve (|Chaboyer et al. 1998j G ratton 
et al. 1997). Integrating over 17m and J7a, the best fit value 
of the age in Hubble-time units is Htfo = 0.93^ ^ or equiv- 
alently f () = 14.5+[;g (0.63/^) Gyr. The age would be some- 
what larger in a flat universe: Htf^ 1 = 0.96+QQ7 or equivalently 



,flat . 



14.9^;t(0.63//i) Gyr. 



Third, even if the universe is not flat, the confidence regions 
of Figure [7] suggest that the cosmological constant is a signifi- 
cant constituent of the energy density of the universe. The best- 
fit model (the center of the shaded contours) indicates that the 
energy density in the cosmological constant is ^0.5 more than 
that in the form of mass energy density. All of the alternative 
fits listed in Table [}] indicate a positive cosmological constant 
with confidence levels of order 99%, even with the systematic 
uncertainty included in the fit or with a clumped-matter metric. 

Given the potentially revolutionary nature of this third con- 
clusion, it is important to reexamine the evidence carefully to 
find possible loopholes. None of the identified sources of statis- 
tical and systematic uncertainty described in the previous sec- 
tions could account for the data in a A = universe. If the uni- 
verse does in fact have zero cosmological constant, then some 
additional physical effect or "conspiracy" of statistical effects 
must be operative — and must make the high-redshift super- 
novae appear almost 0.15 mag (~15% in flux) fainter than the 
low-redshift supernovae. At this stage in the study of SNe la, 
we consider this unlikely but not impossible. For example, as 
mentioned above, some carefully constructed smooth distribu- 
tion of large-grain-size grey dust that evolves similarly for el- 
liptical and spiral galaxies could evade our current tests. Also, 
the full dataset of well-studied SNe la is still relatively small, 
particularly at low redshifts, and we would like to see a more 
extensive study of SNe la in many different host-galaxy envi- 
ronments before we consider all plausible loopholes (including 
those listed in Table ^B) to be closed. 

Many of these residual concerns about the measurement 
can be addressed with new studies of low-redshift supernovae. 
Larger samples of well-studied low-redshift supernovae will 
permit detailed analyses of statistically significant SN la sub- 
samples in differing host environments. For example, the 
width-luminosity relation can be checked and compared for su- 



pernovae in elliptical host galaxies, in the cores of spiral galax- 
ies, and in the outskirts of spiral galaxies. This comparison 
can mimic the effects of finding high-redshift supernovae with 
a range of progenitor ages, metallicities, etc. So far, the re- 
sults of such studies with small statistics has not shown any dif- 
ference in width-luminosity relation for this range of environ- 
ments. These empirical tests of the SNe la can also be comple- 
mented by better theoretical models. As the datasets improve, 
we can expect to learn more about the physics of SN la ex- 
plosions and their dependence on the progenitor environment, 
strengthening the confidence in the empirical calibrations. Fi- 
nally, new well-controlled, digital searches for SNe la at low 
redshift will also be able to further reduce the uncertainties due 
to systematics such as Malmquist bias. 

6.1. Comparison with Previous Results 

A comparison with the first supernova measurement of the 
cosmological parameters in P97 highlights an important aspect 
of the current measurement. As discussed in Section 3, the P97 
measurement was strongly skewed by SN 1994H, one of the 
two supernovae that are clear statistical outliers from the cur- 
rent 42-supernova distribution. If SN 1994H had not been in- 
cluded in the P97 sample, then the cosmological measurements 
would have agreed within the ler error bars with the current 
result. (The small changes in the Zf-corrections discussed in 
Section 2 are not a significant factor in arriving at this agree- 
ment.) With the small P97 sample size of seven supernovae 
(only five of which were used in the P97 width-corrected anal- 
ysis), and somewhat larger measurement uncertainties, it was 
not possible to distinguish SN 1994H as the statistical outlier. 
It is only with the much larger current sample size that it is easy 
to distinguish such outliers on a graph such as Figure ||(c). 

The fact that there are any outliers at all raises one caution- 
ary flag for the current measurement. Although neither of the 
current two outliers is a clearly aberrant SN la (one has no 
SN la spectral confirmation and the other has a relatively poor 
constraint on host-galaxy extinction), we are watching care- 
fully for such aberrant events in future low- and high-redshift 
datasets. Ideally, the one-parameter width-luminosity relation- 
ship for SNe la would completely account for every single well- 
studied SN la event. This is not a requirement for a robust mea- 
surement, but any exceptions that are discovered would provide 
an indicator of as-yet undetected parameters within the main 
SN la distribution. 

Our firs t presentation of the co smological parameter mea- 
surement (Perlmutter et al. 1998c), based on 40 of the current 



42 high-redshift supernovae, found the same basic results as the 
current analysis: A flat universe was shown to require a cosmo- 
logical constant, and only a small region of low-mass-density 
parameter space, with all the systematic uncertainty included, 
could allow for A = 0. (Fit M of Figure ||(f) still shows almost 
the same confidence region, with the same analysis approach). 
The current confidence region of Figure [7] has changed very 
little from t he corresponding confidence region of Perlmutter 
et al. (1998c), but since most of the uncertainties in the low- 
redshift dataset are now included in the statistical error, the re- 
maining systematic error is now a small part of the error budget. 

The more recent a nalysis of 16 high-redshift supernovae 
by Riess et al. (1998) also show a very similar Om-^a con- 
fidence region. The best fits for mass density in a flat- 



universe are QJf = 0.28 ± 0.10 or ft™ = 0.16 ± 0.09 for the 
two alternative analyses of their 9 independent, well-observed, 
spectroscopically-confirmed supernovae. The best fits for the 
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age of the universe for these analy ses are Hptp = 0.90^ ^ and 
Hptp = 0.98^005. To first order, the Riess et al. result provides 
an important independent cross-check for all three conclusions 
discussed above, since it was based on a separate high-redshif t 
supernova search and analysis chain (see Schmidt et al. 1998 ). 
One caveat, however, is that their 51m-^a confidence-region re- 
sult cannot be directly compared to ours to check for indepen- 
dent consistency, because the low-redshift-supernova datasets 
are not independent: a large fraction of these supernovae with 
the highest weight in both analyses are from the Calan/Tololo 
Supernova Survey (which provided many well-measured super- 
novae that were far enough into the Hubble flow so that their pe- 
culiar velocities added negligible redshift-uncertainty). More- 
over, two o f the 16 high-redshift supernovae included in the 



Riess et al. confidence-region analyses were from our sample 



of 42 Supernova Cosmology Project supernovae; Riess et al. 
included them with an alternative analysis technique applied 
to a subset of our photometry results. (In particular, their re- 
sult uses the highest-redshift supernova from our 42-supernova 
sample, which has strong weight in our analysis due to the ex- 
cellent Hubble Space Telescope photometry.) Finally, although 
the analysis techniques are mostly independent , the K correc- 
tions are based on the same Mugent et al. (1998| ) approach dis- 
cussed above. 



6.2. Comparison with Complementary Constraints 
on Qm and J7a 

Significant progress is being made in the measurement of the 
cosmological parameters using complementary techniques that 
are sensitive to different linear combinations of f^M and 
and have different potential systematics or model dependencies. 
Dynamical methods, for example, are particularly sensitive to 
flyi, since J7a affects dynamics only weakly. Since there is 
evidence that dynamical estimates of £Im depend on scale, the 
most appropriate measures to compare with our result are those 
obtained on large scales. From the abundance — indeed the 
mere existence — of rich clusters at high redshift, B ahcall & Fan 
(1998) find n M = 0.2!;; f (95% confidence). The CNOC collab- 



oration ( |Carlberg et al. 19% , 1998) apply evolution-corrected 
mass-to-light ratios determined from virial mass estimates of 
X-ray clusters to the luminosity density of the universe and find 
fl M = 0.17 ±0.07 for il A = (~90% confidence), with small 
chan ges in these results for different values of Oa (cf. C arlberg 
1997 ). Detailed studies of the peculiar velocities of galaxies 
(e.g., [Willick et alj |1997|; |Willick & Strauss ||1998| ; |Riess et al " 



1997b 



but see Sigad et al. 1998) are now giving estimates of 



P=n^/b 1RAS « 0.45 ±0. 1 1 (95% confidence) 1 , where b is the 
ratio of density contrast in galaxies compared to that in all mat- 
ter. Under the simplest assumption of no large-scale biasing for 
IRAS galaxies, b=l, these results give M « 0.26±0.1 1 (95% 
confidence), in agreement with the other dynamical estimates — 
and with our supernova results for a flat cosmology. 

A form of the angular-size distance cosmological test has 



been developed in a series of papers (cf. 3uerra & Daly 1998 



and references therein) and implemented for a sample of four- 
teen radio galaxies by Daly, Guerra, & Wan (1998). The 
method uses the mean observed separation of the radio lobes 
compared to a canonical maximum lobe size — calculated from 
the inferred magnetic field strength, lobe propagation velocity, 
and lobe width — as a calibrated standard ruler. The confidence 



regio n in the £1m-^a plane shown in Daly, Guerra, & Wan 
(1998) is in broad agreement with the SN la results we report; 
they find S1m = 0.2+Q2 (68% confidence) for a flat cosmology. 

QSO gravitational lensing statistics are dependent on both 
volume and relative distances, and thus are more sensitive to 
Using gravitational lensing statistics, ECochanek (1996) 



finds Q\ < 0.66 (at 95% confidence for £1m + £1a = 1), and 
£1 M > 0. 15. ^alco, Kochanek, & Munoz (1998| ) obtained further 
information on the redshift distribution of radio sources which 
allows calculation of the absolute lensing probability for both 
optical and radio lenses. Formally their 90% confidence lev- 
els in the 17m-^a plane have no overlap with those w e report 
here. However, as Falco, Kochanek, & Munoz (1998) discuss, 



these results do depend on the choi ce of galaxy sub-type l umi- 
nosity functions in the lens models. Chiba & Yoshii (1998 ) em- 
phasized this point, reporting an analysis with E/S0 luminosity 
functions that yielded a best-fit mass density in a flat cosmology 
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in agreement with our SN la results. 



Several papers have emphasized that upcoming balloon and 
satellite studies of the Cosmic Background Radiation (CBR) 
should provide a good measurement of the sum of the energy 
densities, SIm + ^a, and thus provide al most orthogo nal infor- 
mation to t he supernova measurements dWhite 1998j T egmark 
et al. 1998). In particular, the position of the first acoustic 
peak in the CBR power spectrum is sensitive to this combi- 
nation of the cosmological parameters. The current results, 
while not conclusive, are already somewhat inconsistent with 
over-closed (f^M + ^A » 1) cosmologies and "near-empty" 
(£1m + ^a <, 0.4) cosmologies, and may excl ude the upper right 
and lower left regions of Figure [7] (see, e.g., Lineweaver 1998 
Efstathiou et al. 1998b. 



6.3. Cosmological Implications 

If, in fact, the universe has a dominant energy contribution 
from a cosmological constant, there are two coincidences that 
must be addressed in future cosmological theories. First, a cos- 
mological constant in the range shown in Figure [7] corresponds 
to a very small energy density relative to the vacuum-energy- 
density scale of particle- physics energy zero-points (see C ar- 
roll, Press, & Turner 1992, for a discussion of this point). Pre- 
viously, this had been seen as an argument for a zero cosmolog- 
ical constant, since presumably some symmetry of the particle- 
physics model is causing cancelations of this vacuum energy 
density. Now, it would be necessary to explain how this value 
comes to be so small, yet non-zero. 

Second, there is the coincidence that the cosmological con- 
stant value is comparable to the current mass-energy density. 
As the universe expands, the matter energy density falls as the 
third power of the scale, while the cosmological constant re- 
mains unchanged. One therefore would require initial con- 
ditions in which the ratio of densities is a special, infinitesi- 
mal value of order 10~ 100 in order for the two densities to co- 
incide today. (The cross-over between mass-dominated and 
A-dominated energy density occurred at z ~ 0.37, for a flat 
£1m ~ 0.28 universe, whereas the cross-over between decelera- 
tion and acceleration occurred when (1 +z) 3 J1m/2 = ft\, that is 
at z ~ 0.73. This was approximately when SN 1997G exploded, 
over 6 billion years ago.) 

It has been suggested that these cosmological coincidences 
could be explained if the magnitude-redshift relation we find 
for SNe la is due not to a cosmological constant, but rather to 



This is an enacjEeighled mpan of Willick et al. (1997) and Riess et al. (1997b), with optical results converted to equivalent IRAS results using &Opt/^lRAS : 
1.20±0.05 frombliver etal. (199f )■ 
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a different, previously unknown physical entit y that contributes 



to the universe's total energy density (see, e.g., Steinhardt 1996 



Turner & White 1997; Caldwell, Dave, & Steinhardt 1998) 



Such an entity can lead to a different expansion history than 
the cosmological constant does, because it can have a different 
relation ("equation of state") between its density p and pressure 
p than that of the cosmological constant, pa/ PA = — 1. We can 
obtain constraints on this equation-of-state ratio, w = p/p, and 
check for consistency with alternative theories (including the 
cosmological constant with w = -1) by f itting the alternative 
expansion histories to data; White (1998) has discussed such 
constraints from earlier supernovae and CBR results. In Fig- 
ure [l0| we update these constraints for our current supernova 
dataset, for the simplest case of a flat u niverse and an equation 
of state that does not vary in time (cf. 3arnavich et al. 1998b, 
for comparison with their high-redshift supernova dataset, and 
Aldering et al. 1998 for time-varying equations of state fit to 
our dataset). In this simple case, a cosmological-constant equa- 
tion of state can fit our data if the mass density is in the range 
0.2 < £Im < 0.4. However, all the cosmological models shown 
in Figure [l(^ still require that the initial conditions for the new 
energy density be tu ned with extreme precision to reach their 
current-day values. Zlatev, Wang, & Steinhardt (1998) have 
shown that some time-varying-w theories naturally channel the 
new energy density term to "track" the matter term, as the uni- 
verse expands, leading — without coincidences — to values of an 
effective vacuum energy density today that are comparable to 
the mass energy density. These models require w J> -0.8 at all 
times up to the present, for f^M > 0.2. The supernova dataset 
presented here and future complementary datasets will allow us 
to explore these possibilities. 
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APPENDIX 



EXTINCTION CORRECTION USING A BAYESIAN PRIOR 



Bayes Theorem provides a means of estimating the a posteriori probability distribution, P(A \A m ), of a variable A given a measure- 
ment of its value, A m , along with a priori information, P(A), about what values are likely: 

mAm) ~ JP(A m \A)P(A)dA (A1) 

In practice P(A) often is not well known, but must be estimated from sketchy, and possibly biased, data. For our purposes here 
we wish to distinguish between the true probability distribution, P(A), and its estimated or assumed distribution, often called the 
Bayesian prior, which we denote as "P(A). Riess, Press, & Kirshner (1996; RPK) present a Bayesian method of correcting SNe la for 
host galaxy extinction. For V{A) they assume a one-sided Gaussian function of extinction, Q{A), with dispersion <jg = \ magnitude: 




for A >0 

V(A) = G(A) ={ V na S (A2) 
forA<0 

which reflects the fact that dust can only redden and dim the light from a supernova. The probability distribution of the measured 
extinction, A,„, is an ordinary Gaussian with dispersion a m , i.e., the measurement uncertainty. RPK choose the most probable value 
of P(A\A m ) as their best estimate of the extinction for each supernova: 
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A m <Tn 



for A,„ > 



Ag = mode(P(A \A m )) ={ a S + a m (A3) 

forA m <0 

Although this method provides the best estimate of the extinction correction for an individual supernova provided V(A) = P(A), 
once measurement uncertainties are considered its application to an ensemble of SNe la can result in a biased estimate of the ensemble 
average extinction whether or not V(A) = P(A). An extreme case which illustrates this point is where the true extinction is zero for 
all supernovae, i.e., P(A) is a delta function at zero. In this case, a measured value of E(B-V) < (too blue) results in an extinction 
estimate of Ag = 0, while a measured value with E(B-V) > results in an extinction estimate Ag > 0. The ensemble mean of these 
extinction estimates will be 

(Ag) = -^(-j^t), CM) 

rather than as it should be. (This result is changed only slightly if the smaller uncertainties assigned to the least extincted SNe la 
are incorporated into a weighted average.) 

The amount of this bias is dependent on the size of the extinction-measurement uncertainties, a m = 1Zb<Je(b-v)- For our sample of 
high-redshift supernovae, typical values of this uncertainty are a m ~ 0.5, while for the low-redshift supernovae, a m ~ 0.07. Thus, if 
the true extinction distribution is a delta-function at A = 0, while the one-sided prior, G(A), of Equation |A2| is used, the bias in (Ag) 
is about 0.13 mag in the sense that the high-redshift supernovae would be overcorrected for extinction. Clearly, the exact amount 
of bias depends on the details of the dataset (e.g., color uncertainty, relative weighting), the true distribution P(A), and the choice 
of prior V(A). This is a worst-case estimate, since we believe that the true extinction distribution is more likely to have some tail 
of events with extinction. Indeed, numerical calculations using a one-sided Gaussian for the true distribution, P(A), show that the 
amount of bias decreases as the Gaussian width increases away from a delta function, crosses zero when P(A) is still much narrower 
than V(A), and then increases with opposite sign. One might use the mean of P(A\A m ) instead of the mode in equation A-3, since the 
bias then vanishes if V(A) = P(A), however this mean-calculated bias is even more sensitive to V{A) ^P{A) than the mode-calculated 
bias. 

We have only used conservative priors (which are somewhat broader than the true distribution, as discussed in Section 4.1), 
however it is instructive to consider the bias that results for a less conservative choice of prior. For example, an extinction distribution 
with only half of the supernovae distributed in a one-s ided Gaussian and half in a delta function at zero extinction is closer to the 



simulations given by Hatano, Branch, & Deaton (1997). The presence of the delta-function component in this less conservative prior 
assigns zero extinction to the vast majority of supernovae, and thus cannot produce a bias even with different uncertainties at low 
and high redshift. This will lower the overall bias, but it will also assign zero extinction to many more supernovae than assumed in 
the prior, in typical cases in which the measurement uncertainty is not significantly smaller than the true extinction distribution. A 
restrictive prior, i.e. one which is actually narrower than the true distribution, can even lead to a bias in the opposite direction from a 
conservative prior. 

It is clear from Bayes Theorem itself that the correct procedure for determining the maximum-likelihood extinction, A, of an 
ensemble of supernovae is to first calculate the a posteriori probability distribution for the ensemble: 

P(A\{A U-J&R. (A5 ) 
( K A "" }) - fP(A)J2P(A mi \A)dA (A5) 

and then take the most probable value of P(A|{A mi }) for A. For the above example of no reddening, this returns the correct value of 
A = 0. 

In fitting the cosmological parameters generally one is not quite as interested in the ensemble extinction as in the combined 
impact of individual extinctions. In this case P(A\{A mi }) must be combined with other sources of uncertainty for each supernova in 
a maximum-likelihood fit, or the use of a Bayesian prior must be abandoned. In the former case a % 2 fit is no longer appropriate 
since the individual P(A\{A mj }y& are strongly non-Gaussian. Use of a Gaussian uncertainty for Ag based on the second-moment of 
P(A|{A mi }) may introduce additional biases. 
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0.19 




1997Q 
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0.03 


0.09 
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22.57 


0.18 




1997R 
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0.001 


23.28 


0.05 


0.11 


-0.66 


23.83 


23.83 


0.23 




1997S 


0.612 


0.001 


23.03 


0.05 


0.11 


-0.67 


23.59 


23.69 


0.21 




1997ac 


0.320 


0.010 


21.38 


0.03 


0.09 


-0.55 


21.83 


21.86 


0.18 




1997af 


0.579 


0.001 


22.96 


0.07 


0.09 


-0.68 


23.54 


23.48 


0.22 




1997ai 


0.450 


0.010 


22.25 


0.05 


0.14 


-0.71 


22.81 


22.83 


0.30 


H 


1997aj 


0.581 


0.001 


22.55 


0.06 


0.11 


-0.68 


23.12 


23.09 


0.22 




1997am 


0.416 


0.001 


21.97 


0.03 


0.11 


-0.67 


22.52 


22.57 


0.20 




1997ap 


0.830 


0.010 


23.20 


0.07 


0.13 


-1.23 


24.30 


24.32 


0.22 


H 



Col 1 . IAU Name assigned to SCP supernova. 

Col 2. Geocentric redshift of supernova or host galaxy. 

Col 3. Redshift uncertainty. 

Col 4. Peak magnitude from lightcurve fit in observed band corresponding to restframe B-band (i.e. m^ edk = m^ edk or m^ e ' ik ). 
Col 5. Uncertainty in fit peak magnitude. 

Col 6. Galactic extinction in observed band corresponding to restframe B-band (i.e., A x = A R or A/); an uncertainty of 10% is 
assumed. 

Col 7. Representative ^-correction (at peak) from observed band to B-band (i.e., K B x = K B r or K B i); an uncertainty of 2% is 
assumed. 

Col 8. B-band peak magnitude. 

Col 9. Stretch-luminosity corrected effective B-band peak magnitude: m| ffectlve = m^ edk + a(s - 1)-K BX -A X . 

Col 10. Total uncertainty in corrected B-band peak magnitude. This includes uncertainties due to the intrinsic luminosity dispersion 
of SNe Ia of 0.17 mag, 10% of the Galactic extinction correction, 0.01 mag for TT-corrections, 300 km s" 1 to account for peculiar 
velocities, in addition to propagated uncertainties from the lightcurve fits. 
Col 1 1 . Fits from which given supernova was excluded. 
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Table 2 
Calan Tololo SNe Ia Data 





z 




peak 
m obs 


peak 
a obs 


A 


X-BB 


m P eak 
m B 


m B 


(7 m corr 


in otes 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(ii) 


1990O 


0.030 


0.002 


16.62 


0.03 


0.39 


-0.00 


16.23 


16.26 


0.20 




1990af 


0.050 


0.002 


17.92 


0.01 


0.16 


+0.01 


17.75 


17.63 


0.18 




1992P 


0.026 


0.002 


16.13 


0.03 


0.12 


-0.01 


16.02 


16.08 


0.24 




1992ae 


0.075 


0.002 


18.61 


0.12 


0.15 


+0.03 


18.43 


18.43 


0.20 




1992ag 


0.026 


0.002 


16.59 


0.04 


0.38 


-0.01 


16.22 


16.28 


0.20 




1992al 


0.014 


0.002 


14.60 


0.01 


0.13 


-0.01 


14.48 


14.47 


0.23 




1992aq 


0.101 


0.002 


19.29 


0.12 


0.05 


+0.05 


19.19 


19.16 


0.23 




1992bc 


0.020 


0.002 


15.20 


0.01 


0.07 


-0.01 


15.13 


15.18 


0.20 




1992bg 


0.036 


0.002 


17.41 


0.07 


0.77 


+0.00 


16.63 


16.66 


0.21 




1992bh 


0.045 


0.002 


17.67 


0.04 


0.10 


+0.01 


17.56 


17.61 


0.19 




1992bl 


0.043 


0.002 


17.31 


0.07 


0.04 


+0.01 


17.26 


17.19 


0.18 




1992bo 


0.018 


0.002 


15.85 


0.02 


0.11 


-0.01 


15.75 


15.61 


0.21 


B-L 


1992bp 


0.079 


0.002 


18.55 


0.02 


0.21 


+0.04 


18.30 


18.27 


0.18 




1992br 


0.088 


0.002 


19.71 


0.07 


0.12 


+0.04 


19.54 


19.28 


0.18 


B-L 


1992bs 


0.063 


0.002 


18.36 


0.05 


0.09 


+0.03 


18.24 


18.24 


0.18 




1993B 


0.071 


0.002 


18.68 


0.08 


0.31 


+0.03 


18.34 


18.33 


0.20 




19930 


0.052 


0.002 


17.83 


0.01 


0.25 


+0.01 


17.58 


17.54 


0.18 




1993ag 


0.050 


0.002 


18.29 


0.02 


0.56 


+0.01 


17.71 


17.69 


0.20 





Col 1 . IAU Name assigned to Calan Tololo supernova. 

Col 2. Redshift of supernova or host galaxy in Local Group restframe. 

Col 3. Redshift uncertainty. 

Col 4. Peak magnitude from lightcurve fit, in observed B-band. Note that the template lightcurve used in the fit is not identical to the 
template lightcurve used by Hamuy et al. so the best-fit peak magnitude may differ slightly. 
Col 5. Uncertainty in fit peak magnitude. 

Col 6. Galactic extinction in observed B-band; an uncertainty of 10% is assumed. 

Col 7. Representative Zf-correction from observed B-band to restframe B-band; an uncertainty of 2% is assumed. 
Col 8. B-band peak magnitude. 

Col 9. Stretch-luminosity corrected B-band peak magnitude. 

Col 10. Total uncertainty in corrected B-band peak magnitude. This includes uncertainties due to the intrinsic luminosity dispersion 
of SNe Ia of 0.17 mag, 10% of the Galactic extinction correction, 0.01 mag for Zf-corrections, 300 km s" 1 to account for peculiar 
velocities, in addition to propagated uncertainties from the lightcurve fits. 
Col 1 1 . Fits from which given supernova was excluded. 
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Table 3 
Fit Results 



Fit N 


x 2 


DOF 




P(V A > 0) 


Best Fit 


Fit Description 


Inclusive Fits 
A 60 


98 


56 




0.9984 


0.83,1.42 


All SNe 


B 56 


60 


52 


0-26^08 


0.9992 


0.85,1.54 


Fit A, but excluding 2 residual 
outliers and 2 stretch outliers 


Primary Fit 
C 54 


56 


50 


28 +0 09 


9979 

Ui77 / 7 


7^ 1 M 


T^it T4 hut iilcA pvr lnHincr 

2 likely reddened 


Comparison Analysis Techniques 
D ' 54 53 51 


25+ 010 


0.9972 


0.76,1.48 


No stretch correction 3 


E 53 


62 


49 


29+ 012 


0.9894 


0.35,0.76 


Bayesian one-sided extinction corrected 13 


Effect of Reddest SNe 
F 51 59 


47 


26 +0 09 


0.9991 


0.85,1.54 


Fit B SNe with colors measured 


G 49 


56 


45 


90+0.09 

U - Z °-0.08 


0.9974 


0.73,1.32 


Fit C SNe with colors measured 


H 40 


33 


36 


n+o- 11 

u ' J1 -0.09 


0.9857 


0.16,0.50 


Fit G, but excluding 7 next reddest 
and 2 next faintest high-redshift SNe 


Systematic Uncertainty Limits 
I 54 56 50 


24+ 09 

U -^M).08 


0.9994 


0.80,1.52 


Fit C with +0.03 mag systematic offset 


J 54 


57 


50 


^^+o io 

U,JJ -0.09 


0.9912 


0.72,1.20 


Fit C with -0.04 mag systematic offset 


Clumped Matter Metrics 
K 54 57 


50 


^S+ 012 

u - JJ -o.io 


0.9984 


2.90,2.64 


Empty beam metric 


L 54 


56 


50 


34+°- 10 

u --^-0.09 


0.9974 


0.94,1.46 


Partially filled beam metric 



a 0.24 mag intrinsic SNe la luminosity dispersion assumed. 

b Bayesian method of Riess et al. (1 996) with conservative prior (see text and Appendix A) and 0.10 mag intrinsic SNe la luminosity 
dispersion. 

c Assumes additional Bayesian prior of Om < 3, ft a < 3. 
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Table 4 

Summary of Uncertainties and Cross-Checks 



Uncertainty 3 on 
(f>flf,ftf) = (0.28, 0.72) 



(A) Calculated Identified Uncertainties 

Statistical Uncertainties (see §5) 

High-redshift SNe 
Low-redshift SNe 
Total 



<7 Qlial = 

M,A 

0.05 
0.065 
0.085 



Systematic Uncertainties from Identified Entities/Processes 

Dust that reddens (see §4.1.2) <0.03 

i.e., 7^ = 0.5) < 2TZ B (z = 0) 
Malmquist bias difference (see §4.2) <0.04 
K-correction uncertainty (see §2, §3) <0.025 

including zero-points 
Evolution of average (see §4.4) <0.01 

SN la progenitor metallicity 

affecting rest-frame B spectral features 
Total 0.05 



(B) Uncertainties Not Calculated 

Proposed/Theoretical Sources of Systematic Uncertainties 
Evolving grey dust (see §4.4, §4. 1 .3) 
i.e., 7^(2 = 0.5) > 2TZ B (z = 0) 
Clumpy grey dust (see §5) 

SN la evolution effects (see §4.4) b 

Shifting distribution of progenitor mass, 
metallicity, C/O ratio 



Bounds and Tests (see text) 

Test with >3-filter color measurements. 

Would increase SN mag residual dispersion with z- 

Test that spectra match on appropriate date for all z- 
Compare low- and high-redshift lightcurve rise-times, 
and lightcurve timescales before and after maximum. 
Test width-luminosity relation for low-redshift SNe 
across wide range of environments. 
Compare low- and high-redshift subsets from ellipti- 
cals/spirals, cores/outskirts, etc. 



(C) Cross Checks 

Sensitivity to (see §4.5) 

Width-luminosity relation 
Non-SN la contamination 
Galactic extinction model 



"m,a 

<0.03 
<0.05 
<0.04 



Gravitational lensing (see §4.3) 
by clumped mass 



<0.06 



a For the redshift distribution of supernovae in this work, uncertainties in ^t^\ correspond approximately to a factor of 1.3 times 
uncertainties in the relative supernova magnitudes. For ease of comparisons, this table does not distinguish the small differences 
between the positive and negative error bars; see Table |3| for these. 

b The comparison of low- and high-redshift lightcurve rise-times discussed in Section 4.4 theoretically limits evolutionary changes 
in the zero-point of the lightcurve width-luminosity relation to less than ~0.1 mag, i.e. A^nm <J 0.13. 
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(O m O a ) — 




0.02 0.05 0.1 0.2 0.5 1.0 



redshift Z 

FIG. 1. — Hubble diagram for 42 high-redshift Type la supernovae from the Supernova Cosmology Project, and 18 low-redshift Type la supernovae from the 
Calan/Tololo Supernova Survey, after correcting both sets for the SN la lightcurve width-luminosity relation. The inner error bars show the uncertainty due to 
measurement errors, while the outer error bars show the total uncertainty when the intrinsic luminosity dispersion, 0.17 mag, of lightcurve- width-corrected Type la 
supernovae is added in quadrature. The unfilled circles indicate supernovae not included in Fit C. The horizontal error bars represent the assigned peculiar velocity 
uncertainty of 300 km s ■ The solid curves are the theoretical »tj ctlVB fe) for a range of cosmological models with zero cosmological constant: (Hm.^a) = (0,0) 
on top, (1,0) in middle and (2,0) on bottom. The dashed curves are for a range of flat cosmological models: (Hm^a) = (0, 1) on top, (0.5,0.5) second from top, 
(1,0) third from top, and (1.5,-0.5) on bottom. 
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Fig. 2.— 
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(a) Hubble diagram for 42 high-redshift Type la supernovae from the Supernova Cosmology Project, and 18 low-redshift Type 
la supernovae from the Calan/Tololo Supernova Survey, plotted on a linear redshift scale to display details at high redshift. The 
symbols and curves are as in Figure 1 . (b) The magnitude residuals from the best-fit flat cosmology for the Fit C supernova subset, 
(£! m ,Oa) = (0.28,0.72). The dashed curves are for a range of flat cosmological models: (Om^a) = (0, 1) on top, (0.5,0.5) third 
from bottom, (0.75,0.25) second from bottom, and (1,0) is the solid curve on bottom. The middle solid curve is for (f^M, ^a) = (0,0). 
Note that this plot is practically identical to the magnitude residual plot for the best-fit unconstrained cosmology of Fit C, with 
(£! m ,^a) = (0.73, 1.32). (c) The uncertainty-normalized residuals from the best-fit flat cosmology for the Fit C supernova subset, 
(O m ,O a ) = (0.28,0.72). 
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FIG. 3. — The distribution of restframe B-band magnitude residuals from the best-fit flat cosmology for the Fit C supernova subset, for (a) 18 Calan/Tololo 
supernovae, at redshifts z < 0.1 and (b) 42 supernovae from the Supernova Cosmology Project, at redshifts between 0.18 and 0.83. The darker shading indicates 
those residuals with uncertainties less than 0.35 mag, unshaded boxes indicate uncertainties greater than 0.35 mag, and dashed boxes indicate the supernovae that are 
excluded from Fit C. The curves show the expected magnitude residual distributions if they are drawn from normal distributions given the measurement uncertainties 
and 0.17 mag of intrinsic SN la dispersion. The low-redshift expected distribution matches a Gaussian with a = 0.20 mag (with error on the mean of 0.05 mag), 
while the high-redshift expected distribution matches a Gaussian with a = 0.22 mag (with error on the mean of 0.04 mag). 
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FIG. 4. — The distribution of lightcurve widths for (a) 18 Calan/Tololo supernovae, at redshifts z < 0.1 and (b) 42 supemovae from the Supernova Cosmology 
Project, at redshifts between 0.18 and 0.83. The lightcurve widths are characterized by the "stretch factor" t, thai stretches prxoattads thp timp axis of a template 



SN la lightcurve to best fit the observed lightcurve for each supernova (see Perlmutter et al. 1995a, 1997e 
been trme-riilateri by a factor 1+7 before fitting to the ohserveri lightcnrves to account for the msmologica 



lengthe 



Kim et al. 1998 ; Kjoldhaber et al. 
g of the" 



199: 



■ supernova timescnle 



The template has 
XGoldhaber et al. 



1995; Leibundgut et al. 1996i ). the shading indicates those measurements ot * with uncertainties less than (J.I, andthe dashed lines indicate the 'two supernovae 
that are removed trom the tits after Fit A. These two excluded supernovae are the most significant deviations from s = I (the highest-stretch supernova in panel 
(b) has an uncertainty of ±0.23 and hence is not a significant outlier from s = 1); the remaining low- and high-redshift distributions have almost exactly the same 
error- weighted means: (j')Hamuy = 0.99 ± 0.01 and (s)scp = 1. 00 ± 0.01. 
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Comparison of best-fit confidence regions in the Om-^a plane. Each subpanel shows the result of Fit C (shaded regions) compared 
with fits to different subsets of supernovae, or variant analyses for the same subset of supernovae, to test the robustness of the Fit C 
result. Unless otherwise indicated, the 68% and 90% confidence regions in the Om-^a plane are shown, after integrating the four- 
dimensional fits over the other two variables, A4 B and a. The "no-big-bang" upper-left shaded region, the fiat-universe diagonal 
line, and the infinite expansion line are shown as in Figure [7] for ease of comparison. The subpanels show: (a) Fit A of all 60 
supernovae; and Fit B of 56 supernovae, excluding the two outliers from the lightcurve width distribution and the two remaining 
statistical outliers. Fit C further excludes the two likely reddened high-redshift supernovae. (b) Fit D of the same 54-supernova 
subset as in Fit C, but with no correction for the lightcurve width-luminosity relation, (c) Fit H of the subset of supernovae with color 
measurements, after excluding the reddest 25% (9 high-redshift supernovae) and the two faint high-redshift supernovae with large 
uncertainties in their color measurements. The close match to the confidence regions of Fit C indicates that any extinction of these 
supernovae is quite small, and not significant in the fits of the cosmological parameters, (d) The 68% confidence region for Fit E of 
the 53 supernovae with c olor measurements from the Fit B dataset, but following the Bayesian reddening-correction method of Riess, 
Press, & Kirshner (1996). This method, when used with any reasonably conservative prior (i.e., somewhat broader than the likely 
true extinction distribution; see text), can produce a result that is biased, with an approximate bias direction and worst-case amount 
indicated by the arrows, (e) Fits I and J are identical to Fit C, but with 0.03 or 0.04 magnitudes added or subtracted from each of 
the high-redshift supernova measurements, to account for the full range of identified systematic uncertainty in each direction. Other 
hypothetical sources of systematic uncertainty (see Table f|B) are not included, (f) Fit M is a separate two-parameter (J1m and J7a) fit 
of just the high-redshift supernovae, using the values of Mb and a found from the low redshift supernovae. The dashed confidence 
regions show the approximate range of uncertainty from these two low-redshift-derived parameters, added to the systematic errors of 
Fit J. Future well-observed low-redshift supernovae can constrain this dashed-line range of uncertainty. 
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FIG. 6. — (a) The restframe B—V color excess distribution for 17 of 18 Calan/Tololo supernovae (see text), corrected for Galactic extinction using values from 
Schlegel et. al. 1998. (b) The restframe B—V color excess for the 36 high-redshift supernovae for which restframe B-V colors were measured, also corrected for 
Galactic extinction. The darker shading indicates those E(B-V) measurements with uncertainties less than 0.3 mag, unshaded boxes indicate uncertainties greater 
than 0.3 mag, and the light shading indicates the two supernovae that are likely to be reddened based on their joint probability in color excess and magnitude 
residual from Fit B. The dashed curve shows the expected high-redshift E{B-V) distribution if the low-redshift distribution had the measurement uncertainties of 
the high-redshift supernovae indicated by the dark shading. Note that most of the color-excess dispersion for the high-redshift supernovae is due to the rest-frame 
V-band measurement uncertainties, since the rest-frame B-band uncertainties are usually smaller. 
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Fig. 7. — Best-fit confidence regions in the Qm-^a plane for our primary analysis, Fit C. The 68%, 90%, 95%, and 99% statistical confidence regions in 
the Qm=Qa plane ate shown aftpr integrating_the four-dimensional fit over Mb and a. (The table of this two-dimensional probability distribution is available 
at ittp://www-supernova. lbl.gov/.) See Figure pie) for limits on the small shifts in these contours due to identified systematic uncertainties. Note that the spatial 
curvature of the universe — open, flat, or closed — is not determinative of the future of the universe's expansion, indicated by the near-horizontal solid line. In 
cosmologies above this near-horizontal line the universe will expand forever, while below this line the expansion of the universe will eventually come to a halt and 
recollapse. This line is not quite horizontal because at very high mass density there is a region where the mass density can bring the expansion to a halt before the 
scale of the universe is big enough that the mass density is dilute with respect to the cos molagica] constant pnprgy densi ty The upper-left shaded region, labeled "no 
big bang," represents "bouncing universe" cosmologies with tin hig hang in the past (see [Carroll, Press, & Turner 1992 ). The lower right shaded region corresponds 
to a universe that is younger than the oldest heavy elements ([Schramm 1990|), for any value of Hq > 5U km s _1 Mpc~' ! 
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FIG. 8. — Best-fit 68% and 90% confidence regions in the Qyi-^A plane for cosmological models with small scale clumping of matter (e.g., in the form 
of MACHOs) compared with the Friedman-Robertson- Walker model of Fit C, with smooth small-scale matter distribution. The shaded contours (Fit C) are the 
confidence regions fit to a Friedman-Robertson-Walker magnitude-redshift relation. The extended confidence strips (Fit K) are for a fit of the Fit C supernova set to 
an "empty beam" cosmology, using the "partially filled beam" magnitude-redshift relation with a filling factor r\ — 0, representing an extreme case in which all mass 
is in compact objects. The Fit L un-shaded contours represent a somewhat more realistic partially-filled-beam fit, with clumped matter (r) = 0) only accounting for 
up to Cm = 0.25 of the critical mass density, and any matter beyond that amount smoothly distributed (i.e., r) rising to 0.75 at Qm = 1). 
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planp fnr thp primary analysis Fit C Thp isnrhrnnf 



(Hamiiy et al 'QQ*; Eiess, Eress, & Kkghnet 



arp lahplpfi fnr t hp rase nf H n = 63 km s Mpc , representing a typical value found from studies of SNe la 
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the age labels would be lU7o smaller. The diagonal line labeled accelerating/decelerating is drawn tor qg = Hm/^~ = lJ > an d divides the cosmological 
models with an accelerating or decelerating expansion at the present time. 
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FIG. 10. — Best-fit 68%, 90%, 95%, and 99% confidence regions in the Qm _ w plane for an additional energy density component, Q w , characterized by an 
equation-of- state w = p/p. (If this energy density component is Einstein's cosmological constant, A, then the equation of state is w = p\/ Pa = ""!•) The fit is for 
the supernova subset of our primary analysis, Fit C, constrained to a flat cosmology (SIm +f2 w = 1). The two variables Mb and a are included in the fit, and then 
integrated over to obtain the two-dimensional probability distribution shown. 



